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ABSTRACT 


The  aeroelastic  flutter  and  divergence  behavior  of 
rectangular,  graphite/epoxy /cantilevered  plates  with  vary¬ 
ing  amounts  of  bending- torsion^ stiffness  coupling  is  inves¬ 
tigated  for  incompressible  flow.  A  general  Rayleigh-Ritz 
formulation  is  used  to  calculate  flexibility  influence  coef¬ 
ficients,  static  deflections,  divergence  velocities,  vibra¬ 
tion  frequencies,  and  flutter  velocities.  Flutter  calcula¬ 
tions  are  done  using  the  U-g  method.  Test  plates  were  con¬ 
structed  and  subjected  to  static,  vibration  and  wind  tunnel 
tests.  Wind  tunnel  tests  indicated  static  deflections, 
divergence  instabilities,  bending- torsion  flutter  at  low 
angles  of  attack,  and  stall  flutter  at  high  angles  of  attack. 
Bending  stiffness  and  first  bending  frequencies  showed  good 
agreement  between  theory  and  experiment.  Torsional  stiffness 
and  first  torsion  frequencies  were  not  accurately  predicted. 

”by  the  theory  for  highly  coupled  plates.  Divergence  velocities 
and  reduced  flutter  velocities  showed  reasonable  agreement 
between  theory  and  experiment.  Test  plates  with  varying 
amounts  of  coupling  exhibited  markedly  different  stall  flutter 
characteristics.  „ 
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CHAPTER  I 
INTRODUCTION 


1.1  Background 

Aeroelastic  tailoring  is  broadly  defined  as  the  tech¬ 
nology  to  design  a  lifting  surface  which  exhibits  a  desired 
aeroelastic  response.  Desired  aeroelastic  responses  most 
often  considered  are  maximization  of  flutter  and  divergence 
speeds.  Other  aeroelastic  responses  can  be  important  de¬ 
pending  on  the  application.  They  include  control  reversal 
speed,  camber  changes  as  a  function  of  load  and  speed,  and 
angle  of  attack  changes  as  a  function  of  load.  The  aniso¬ 
tropy  of  advanced  composite  materials  or,  more  specifically, 
the  designers '  ability  to  control  that  anisotropy  by  selec¬ 
tive  lamination,  makes  it  an  attractive  material  for  aero¬ 
elastic  tailoring. 

Aeroelastic  tailoring,  which  exploits  the  advantage 

of  advanced  composites,  has  received  considerable  attention 

in  recent  literature.  N.J.  Krone,  Jr.1-  concluded  that 

forward  swept  wings  without  divergence  or  weight  penalties 

may  be  possible  through  the  use  of  selectively  laminated 

2  3 

advanced  composites.  T. A.  Weisshaar  '  extended,  analyti¬ 
cally,  Krone's  conclusion  to  potentially  practical  wing 
designs.  Weisshaar  concluded  that  the  binding-torsion 
stiffness  coupling  of  anisotropic  advanced  composite 
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materials  was  the  key  to  eliminating  divergence  in  forward 

4 

swept  wings.  V.C.  Sherrer,  T.J.  Hertz,  and  M.H .  Shirk 
conducted  a  series  of  wind  tunnel  tests  using  simple  plate¬ 
like  models  of  a  forward  swept  wing.  These  tests  essentially 
verified  Weisshaar's  conclusion,  and  also  showed  that  existing 
analytical  techniques  (computer  programs)  would  adequately 
predict  the  divergence  dynamic  pressures  for  most  test  con¬ 
ditions  . 

Under  contract  to  the  Air  Force  Flight  Dynamics  Lab¬ 
oratory,  General  Dynamics  Corporation  developed  a  wing 

5 

aeroelastic  synthesis  computer  program  (TSO)  .  The  computer 
program  was  intended  to  be  a  preliminary  design  tool  for 
optimization  of  wings  with  advanced  composite  structural  box 
skins.  To  this  end,  it  used  a  direct  Rayleigh-Ritz  energy 
formulation  to  model  the  structural  deflections.  The  prog¬ 
ram  was  capable  of  optimizing  a  wing  skin  design  for  several 
different  constraints,  simultaneously.  Published  aeroelastic 
tailoring  studies  '  '  using  TSO  have  not  examined,  in  depth, 
the  effect  of  structural  bending-torsion  coupling  on  flutter 
speed  in  the  absence  of  changes  in  other  variables.  Addi¬ 
tionally,  only  very  limited  wind  tunnel  data  are  available 
which  examine  the  effect  of  structural  bending  torsion 
coupling  on  flutter  speed,  or  even  that  verify  the  accuracy 
of  flutter  speeds  generated  by  TSO.  Finally,  the  phenomenon 
known  as  stall  flutter  has  received  only  limited  attention 
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8  9 

as  it  relates  to  metal  lifting  surfaces  '  ,  and  virtually 
none  as  related  to  advanced  composite  lifting  surfaces. 

The  references  cited  in  this  section  by  no  means  com¬ 
prise  a  complete  list  of  work  done  in  this  area.  They  should 
be  viewed  more  as  an  indicator  of  where  emphasis  has  been 
placed.  One  obstacle  to  conducting  research  in  aeroelastic 
tailoring  is  that  much  of  the  work  done  by  the  larger  aero¬ 
space  corporations  is  proprietary  and  consequently  is  never 
reported  in  literature  accessible  to  the  public. 

1.2  Objectives 

This  study  will  attempt  to  ascertain  the  effect  of 
varying  amounts  of  bending- tors ion  coupling  on  both  the 
divergence  and  flutter  speeds  of  an  unswept  lifting  surface 
in  incompressible  flow.  The  lifting  surface  will  be  ideal¬ 
ized  by  a  cantilevered,  rectangular,  flat  plate  constructed 
of  laminated  graphite/epoxy,  and  having  a  half-span  aspect 
ratio  of  4  (8  by  standard  aerodynamic  convention) . 

Specific  objectives  are  to 

(1)  Develop  an  analytical  formulation  to  predict  static 

deflections,  natural  vibration  frequencies,  divergence 
speed,  and  flutter  speed  for  idealized  lifting  surfaces. 
The  formulation  must  be  applicable  to  a  plate  having 
substantial  bending-torsion  coupling,  but  limited  to 
a  mid-plane  symmetric  lamination  arrangement. 


(2)  Evaluate  the  effect  of  bending-torsion  coupling  on 
the  idealized  lifting  surface  by  performing  static 


deflection,  vibration,  and  wind  tunnel  tests  on  sel¬ 
ected  cantilevered  plates  having  varying  amounts  of 
bending-torsion  coupling. 

(3)  Investigate  the  effects  of  bending-torsion  coupling  on 
stall  flutter  velocity.  This  will  be  limited  to  an 
experiment  evaluation. 

(4)  Determine  the  accuracy  of  the  analytical  formulation 
by  comparing  theoretical  and  experimental  results . 

1 . 3  Organization 

Chapter  II  develops  a  general  Rayleigh-Ritz  energy 
formulation  to  model  lateral  deflections  of  a  laminated 
cantilever  plate.  Solutions  for  both  two  term  and  three 
term  deflection  equations  are  derived.  The  two  term  general 
solution  is  than  applied  to  a  static  deflection  problem  to 
calculate  a  2  x  2  matrix  of  flexibility  influence  coef¬ 
ficients,  and  to  a  divergence  problem  to  calculate  the 
divergence  velocity  of  the  idealized  lifting  surface.  Both 
two  and  three  term  solutions  are  applied  to  a  free  vibration 
problem  to  determine  the  lower  two  or  three*  vibration  fre¬ 
quencies.  Finally,  the  two  term  solution  is  applied  to  a 


classical  bending-torsion,  potential  flow,  flutter  problem 
using  the  U-g  method. 

Chapter  III  presents  the  experimental  test  apparatus 
and  procedures  for  the  static  deflection  tests,  free  vibra¬ 
tion  tests,  and  wind  tunnel  tests  (flutter  and  divergence) . 
The  test  results  are  then  discussed,  commenting  on  possible 
sources  of  error. 

Chapter  IV  compares  the  analytical  and  theoretical 
results,  identifies  areas  of  good  and  poor  agreement,  and 
comments  on  sources  of  inaccuracies  in  the  theoretical 
analysis,  as  well  as  ways  to  improve  it. 

Chapter  V  summarizes  the  conclusions  reached  in  the 
previous  chapters  and  makes  recommendations  for  further 
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CHAPTER  II 
THEORY 


2 . 1  Anisotropic  Plate  Flexural  Stiffness 

The  flexural  modulus  (stiffness)  components  of  a  lami¬ 
nated  advanced  composite  (graphite/epoxy  in  this  study)  plate 
are  dependent  on  both  the  fiber  orientation  and  stacking 
sequence  of  the  individual  laminae  (plies>.  To  simplify 
fabrication  of  test  specimens,  only  laminated  plates  with 
a  mid-plane  symmetric  stacking  sequence  are  considered  . 

The  ply  angles  (0)  follow  the  sign  convention  in  Fig.  1. 


Figure  1.  Ply  Angle  Sign  Convention 
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The  in-plane,  on-axis  lamina  modulus  components  (Q^j) 
were  obtained  from  the  orthotropic  engineering  constants 
for  Hercules  AS/3501-6  graphite/epoxy,  from  which  the  test 
specimens  were  to  be  fabricated.  These  engineering  constants 
take  on  different  values  depending  on  whether  they  are  ob¬ 
tained  from  bending  or  stretching  tests.  Engineering  con¬ 
stants  obtained  from  each  type  of  test  appear  in  Appendix  A, 


and  their 

validity 

for  vibration  problems 

is  briefly  dis- 

cussed  in 

Chapter 

IV  (Reference 

10  has  an 

in-depth  treatment 

of  this  subject.) 

The  terms 

are  defined  as11 

Q11 

-  el/(1 

VLTVTL) 

(2-la) 

Q22 

=  Et/(1 

"  VLTVTL) 

(2-lb) 

Q12 

-  Q21  - 

VLTET/(1  "  V 

LTVTL) 

( 2-lc) 

Q66 

=  glt 

( 2-ld) 

where  vTL  =  (ET/EL> VLT 

( 6 ) 

The  off-axis  lamina  modulus  components  (Q^j  )  were 
obtained  by  first  defining  a  set  of  invariants11 


X1  * 


[Q11  +  Q22  +  2Q12]/4 


Z2  = 


[Q11  +  Q22  -  2Q12  +  4Q661/8 


R1  = 


[Qll  -  Q22]/2 


R2  = 


(Q11  +  Q22  -  2Q12  “  4Q661/8 


( 2-2a) 
( 2-2b) 


(2-2c) 


( 2-2d) 
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and  using 

relations 

the  invariants  in  the 

11. 

following 

transformation 

Qu 

s 

I1  +  I2  +  R1cos20  + 

R2cos40 

(2-3a) 

q<|’ 

= 

Ix  +  I  -  R1cos29  + 

R2cos40 

( 2-3b) 

q<«) 

w12 

= 

ri  “  Z2  ~  R2cos49 

(2- 3c) 

= 

i2  -  R2cos40 

(2-3d) 

< 

= 

R^sin20  +  R2sin40 

( 2-3e ) 

q(6) 

^26 

= 

j  R1sin20  -  R2sin40 

(2-3f) 

whe  re  9 

is 

the  ply  angle  (Fig. 

1)  . 

The  flexural  modulus  (D^Jfor  an  n-ply  laminate  with 
arbitrary  ply  angle  orientation  is  obtained  from 


D. .  =  Z  Q*(?k)[(Zv  _  3  )/3]  i;j  =  lf2i6  (2_4) 

JC  =  1  ^ 

where 

0k  =  ply  angle  of  the  (k)  ply 

=  distance  from  the  mid-plane  to  the  upper  sur¬ 
face  of  the  (k)  ^  ply  (positive  above  mid-plane, 
negative  below  mid-plane) 

zk_^  =  distance  from  the  midplane  to  the  lower  surface 
of  the  (k) th  ply 

Flexural  modulus  values  (D^,  D66'  °16^  ^or 
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[+30-/0]  .  [+45-/0]  ,  and  [±45/0]  laminates  appear  in 
Appendix  A.  The  reader  should  observe  that,  using  the 
equations  in  this  section,  the  flexural  moduli  of  a  [ +6 2/° 1 3 
laminate  will  be  the  same  as  those  for  a  [-02/O)s  laminate, 
with  the  exception  of  a  negative  sign  on  Dlg  and  D2g  for  the 
latter  laminate. 

2 . 2  General  Rayleigh-Ri tz  Formulation 

The  direct  Rayleigh-Ri tz  energy  method  is  a  relatively 
simple,  straightforward  approximation  for  the  plate  deflec¬ 
tions,  as  required  for  the  static  deflection, free  vibration, 
divergence,  and  flutter  analyses  in  this  study.  The  Ray¬ 
leigh-Ri  tz  method  also  has  the  advantage  of  showing  the 
effect  of  the  individual  variables  on  the  solution  more 
clearly  than  other  more  accurate  methods,  such  as  finite 
element  analysis.  Whether  or  not  the  results  of  the  Ray- 
leigh-Ritz  analysis  developed  in  this  chapter  correlate 
sufficiently  well  with  experimental  results  is  discussed  in 
Chapter  IV.  To  simplify  computation,  the  "wing"  is  idealized 
by  a  rectangular  cantilevered  flat  plate  with  uniform  thick¬ 
ness.  Further,  to  allow  the  stiffness  properties  of  the 
plate  to  be  more  accurately  depicted,  an  aerodynamic  fairing 
will  not  be  used  over  the  plate. 

The  Rayleigh-Ritz  analysis  begins  by  assuming  a  deflec¬ 
tion  shape  for  the  structure.  If  only  lateral  deflections 


(w)  are  allowed,  the  single  deflection  equation,  written  in 
generalized  coordiates  is: 
n 

w  =  I  Y-i  (x,y)q.  (t)  (2-5) 

i=l  1  1 

where  y^(x,y)  is  the  nondimensional  deflection  or  mode  shape 
of  the  (i)^  mode  (which  must  satisfy  geometric  boundary 
conditions  for  a  cantilevered  plate),  and  q^(t)  is  the  gen- 
eralized  displacement  of  the  (i)^  mode.  The  generalized 
displacement  is  a  function  only  of  time  and  has  units  of 
length.  For  convenience  and  consistency,  the  coordinate 
system  in  Fig.  2  will  be  used  for  all  problems. 

The  deflection  equation  is  simplified  by  assuming 

(1)  the  plate  is  chordwise  rigid  when  undergoing  torsional 
deflections, 

(2)  the  plate  does  not  exhibit  chordwise  bending  for  any 
of  the  problems  to  be  developed, 

(3)  a  single  term  in  the  deflection  equation  for  each  of 
the  desired  deflection  modes  (first  bending  and  first 
torsion;  or  first  bending,  second  bending,  and  first 
torsion)  will  adequately  represent  the  actual  deflec¬ 
tions  of  the  plate.  A  further  requirement  here,  to 
insure  rapid  convergence,  is  that  the  mode  shape  for 
each  term  must  accurately  depict  the  deflection  for 
that  mode.  The  two  term  deflection  equation  can  now 


be  written  as 


i 


Figure  2.  Sign  Conventions 


w 


<j>1(x)4>1(y)q1(t)  +  *a<x)Vy)q2(t) 


(2-6) 


where  <^(x)  and  (y )  are  one  dimensional  mode  shapes  for 

'f'Vi 

the  (i)  mode.  If  the  first  term  represents  the  first 
bending  mode,  then  <J>1(y)  is  simply  a  rigid-body  translation 
($^(y)  =  1)  •  If  the  second  term  represents  the  first  tor¬ 
sion  mode,  then,  by  the  first  assumption,  <f>a(y)  is  a  rigid 
body  rotation  (<j>  (y)  *  y/c) .  Incorporating  these  into  Eq. 

2-6  yields 

w  =  !J)1(x)q1(t)  +  (y/c)  <J>a  (x)  q2  (t)  (2-7) 

We  can  apply  the  same  analogy  to  the  three  term  deflection 
equation,  which  is  written  as 

w  =  <j)1(x)q1(t)  +  <J>2(x)q2(t)  +  (y/c)  (x)  q3  (t)  (2-8) 

In  Fig.  2b,  the  elastic  axis  (EA)  is,  by  definition, 
located  at  the  midchord  of  the  plate  (y  -  0) .  Lateral  def¬ 
lection  of  the  elastic  axis  is  designated  w£  and  rotation 
of  the  plate  about  its  elastic  axis  is  designated  a.  One 
observes  that  wE  and  a  are  not  generalized  coordinates  and, 
in  fact,  q2  in  Eq.  2-7  and  a  do  not  even  have  the  same  units. 
The  units  of  a  are  radians,  while  q2  has  units  of  length. 

To  write  w^  and  a  in  generalized  coordinates  (for  Eq.  2-7) , 
the  following  relations  are  used: 


i 
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WE  =  w  y 2*0  *  ♦  j.Wq^t)  (2-9) 

a  =  ||  =  («ta(x)/c)q2(t)  (2-10) 


The  importance  of  this  transformation  will  become  apparent 
when  aerodynamic  forces  are  applied.  Finally,  from  Fig.  2b, 
Fg  is  a  lateral  force  applied  at  the  elastic  axis  and  Mg  is 
a  moment  about  the  elastic  axis. 

The  strain  energy  (V)  for  a  symmetric  anisotropic  lami- 

.  12 
nate  is 


V  = 


i  / 
z  0 


c 

2 


c 

'7 


lDll<“'xx> 


2D,_w,  w, 

12  'xx  'yy 


D22("'yy)2  + 


4D, ,w,  w,  + 
16  xx  'yy 


4D_,.w,  w, 

26  yy  xy 


+  4D66(w'xy)  1  dydx 


(2-11) 


where  a  comma  denotes  partial  differentiation. 


2.2.1  Two  Term  Deflection  Equation 


w'xx  -  »'i1l  +  I  *1*2  w,xy  = 


w,yy  -  0 


(2-12) 


where:  (  ) '  =  d/dx  and  (  )"  =  d^/dx^ .  Substituting  Eq.  2-12 


into  Eq.  2-11  and  performing  the  chordwise  integration  yields 


e 
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V  = 


qjl  ^  o/(**l)2dx]  + 


-j  CD.  .  x.  0  tvcc  x,  -j 

+  q ,[  f  (<{,n)2dx  +  — ££■  /  (4)  ’  )  2 dx]  (2-1 


2D, 


l2  24 


The  kinetic  energy  expression  for  the  plate  is 


T  = 


\  / 
z  0 


c 

7  2 

/  m  w  dydx 


(2-1 


where  (')  =  d/dt  and  m  is  the  mass/area. 


m  =  PGEfcp 


(2-1] 


where  pG£  is  the  specific  gravity  of  graphite/epoxy  and  t 

P 


the  plate  thickness 


w  »  +  (y/c)4>aq2 


( 2— 1 1 


Performing  the  chordwise  integration,  Eq.  2-14  then  becomes 


T  =  q*  l  X  Q/%idx]+  q2  (  /  ^dx] 


(2-1' 


The  change  in  external  work  (<5We)  can  be  expressed  as 


<5W„ 


c 

J 


f  c!  PzSw  dydx 


(2-11 


where  pz  is  a  distributed  lateral  load  and 


<5W  =  *,59,  + 


(2-1 
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Then,  by  substitution 


6We  =  Sq^  +  6q2Q2 


where 


Q 


1 


l 

f 

0 


P^dydx 


and 


c 

1  V 

_/  p_  *;  <J>  dydx 
cJ  cz  c 

2 


(2-20) 


The  relationship  between  the  work  and  energy  expressions 
is  obtained  through  Lagrange's  equation,  which  is  a  statement 
of  Hamilton's  energy  principle;  the  basic  premise  of  the 
Rayleigh-Ritz  method.  Lagrange's  equation  is13 


d_ 

dt 


3T  )  +  3V 

3<*i  3q± 


3T  =  Qi  (2-21) 

3q± 

(i  «  1,2, . . .N) 


where : 


3 V 
3<3 


-  0 


3V 

3^ 


dc  i  a 

2q  [  -il-  /  («{,••)  2dx]  +  q,[2D.,  /<J>  *  <(,'  dx] 

A  2  0  *  ■L0  0  a 


12  *  -Jl  I2D16  /  «Wxl  +  2^2  /  <♦;> 

+  -w  j  dx! 


cD 


11 


"dx 


From 

with 
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^  =  2q,  [5^  /  dx] 

aqL  1  ^0  1 

i 

4r  (  —  )  =  q,  [  me  /  <}>?  dx] 

dt  3qL  0 


9 T  r  me  -  ,2  -  ■. 

-  2c?2  [  24  J  dx  ] 

3q2  0 


d  ,  3T  .  _  r  me  ,  ,2  ,  , 

dt  a  •  q2  12  nf  dx5 

3q2  0 


Eq.  2-21,  the  general  equations  of  motion  for  the  plate 

a  two  term  assumed  displacement  equation  are 

i  i 

q  [  me  f  $2  dx]  +  q.  [D..C  /  (<j>")2dx] 

i  0  1  L  xx  o 


+  q2  [  2Dlg  f  dx]  =  ( 2-22a) 


V  r!  /  *1 dxl  +  <11  t  2D16  /  *i*i  dxl 

d^ll  ^  2 

+  ^ir/  <♦;>  dx 

4D  1  ? 

+  — ~  /  (♦: )  dx]  =  q2 

^  o  u 


(2-22b) 


To  make  these  equations  more  compact,  nondimensional  expres¬ 
sions  for  the  integrals  are  defined  below.  These  integral 
expressions  may  be  evaluated  exactly,  or  using  any  of  several 
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numerical  integration  techniques,  if  the  integrand  is  too 
complicated. 


X3 

= 

1 

z 

2, 

/ 

0 

( 2-23a) 

X4 

= 

1 

£ 

z 

s 

0 

(2- 2  3b) 

X5 

= 

1 

I 

£ 

f 

0 

<j>^  dx 

(2-23c) 

z 

X6 

r 

z 2 

/ 

0 

dx 

( 2-23d) 

X7 

= 

z3 

£ 

/ 

0 

(0'{)2  dx 

(2-23e) 

X8 

r 

z 

£ 

/ 

0 

(4^)  2  dx 

( 2-23f ) 

X11 

= 

l 3 

£ 

/ 

n 

($;)2  dx 

( 2  —  2  4  g) 

0 


Incorporating  the  expressions  of  Eq.  2-23  into  Eqs .  2-22a 
and  b  yields  the  final  form  of  the  equations  of  motion. 


<*1 

$2 


D11C 


2D 


16 


[mcU4]  +  qL  [-p-  I?]  +  q2[  I6J  =  Q1  (2.24a) 

r  me*  T  i  ,  _  ,  2D16  t  , 

[  TT  X53  +  ql  [  ~TT~  X6  3 


£ - 

.  r  4D66  ,  D11C"  -  .  , 

+  q2  c£  {I8  48DggJl 7  Ill)  3 


EQUATIONS  OP  MOTION 


(2.24b) 


In  Eq.  2-24b,  the  term 
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48FD^“  "11 

represents  the  warping  stiffness,  which  is  inversely  pro¬ 
portional  to  the  square  of  the  plate  aspect  ratio  ( Z/c ) . 
For  the  plates  in  this  study  (AR  =  4),  this  term's  contri¬ 
bution  will  be  small. 


2.2.2  Three  Term  Deflection  Equation 

The  derivation  of  the  equations  of  motion  for  the  three 
term  deflection  equation  is  almost  identical  to  the  two-term 
deflection  equation  derivation.  Therefore,  only  the  final 
expressions  for  V,  T,  5W_,  and  Lagrange's  equation  will  be 
presented.  As  one  would  expect,  there  will  be  three  equa¬ 
tions  of  motion  instead  of  two. 


Deflection  equation: 


w  =  $1(x)q1(t)  +  <j>2(x)q2(t)  +  ^<J>a  (x)  q3  (t)  (2-25) 


Strain  energy: 


'11'*  "  ?  D-i  -i  C  £  7  7 

V  =  [  2~  Qf  (<*,1)  dxlql  +  1  ~T~  /  <4>2}  dxlq2 


Di  I  c  Z  o 

/  .  It  \  , 


2D66  * 


+  1  -IT"  !  (<>Mx  +  —ir  !  U’)2dx]q^ 

0  a  0  a 

z  z 

+  [  Dnd  /  dx]q1q2  +  [  2Dlg  /  ^^dxlq-jq-j 


+  ^  2D16  ^2^adx^q2q3 


(2-26) 
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Kinetic  energy: 


t  =  qf  [  %£  i  ^:dx]  +  q;  [  ^  ;  4>;dx] 


me 


+  I  U  /  4>‘dx]  +  q1q2  [  me  /  ♦1*2dx] 


(2-27) 


Change  in  external  work: 


5We  =  Q16q1  +  Q25q2  +  Q3<Sq3 


(2-28) 


where : 


i  - 

Q1  ■  of  -£/2  d*dx 


°2  * 


f  J  P24>2  dYdx 


Q3  = 


0 


f  c/2  P2  &cx  dydx 
~2 


Before  applying  Lagrange's  equation,  we  eliminate  terms 

with  i  i 

f  1  4>  _ d.x  and  /  4>"  dx 

0  1  Z  0  1 

by  requiring  that  <t^(x)  and  4>2(x)  orthogonal  functions. 
Now,  using  the  same  nondimensional  expressions  for  the  inte¬ 
grals  as  before,  the  equations  of  motion  obtained  from 
Lagrange's  equations  are: 


q1  [  mc2,  I4  ]  +  qx  C 


i7  1  +  q3  t 


xo 


26 


Dnc 


2D 


J6  1 


-  Qi 


(2-29a) 


q2  [  mcl  I12  ]  +  q2  f 


2D, 


I9  ] 


=  Q- 


(2-29b) 


q3  [  —  I5  J  +  q,  I 

J  12  5  1 


2Dlfi  2D., 

— —  Ifi  J  +  q2  [  — ~ 

2D  2  o  2 


l 

4D 


*9  1 


+  q?  t 


cl 


-  (  I3  +  ^-In)]  -  03 

3  48D/._J,2  11  3 


66' 


(2-29c) 


EQUATIONS  OF  MOTION 


where 


=  Jl2  /  <|>"4> '  dx 

2ra 


2 

u  \  ^ 


i1Q  =  z3  ;  (0n)  dx 


12 


1  2 
2  /  <t>2  dx 

X,  o  ^ 


The  functions  used  in  this  study  for  the  mode  shapes 
<j>2,  <j>a)  are  combinations  of  sin,  cos,  sinh,  and  cosh. 
These  functions  are  known  as  ''beam''  modes  because  they 
very  closely  model  the  deflections  of  an  isotropic  Euler 
beam.  The  first  bending  mode  shape  (<^)  and  the  first 
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torsion  mode  shape  (4>a)  were  obtained  from  Ref.  14  and  the 
second  bending  mode  was  obtained  from  Ref.  15.  These  mode 
shapes,  along  with  numerical  values  for  the  nondimensional 
integral  expressions,  appear  in  Appendix  B. 


2.3  Static  Deflection  Problem 


The  static  deflection  problem  is  formulated  as  an 
analytical  model  of  the  experimental  deflection  tests 
described  in  Chapter  III.  For  this  problem,  the  canti¬ 
levered  plate  is  subjected  to  a  concentrated  unit  load  and 
moment,  individually  applied  at  the  3/4  length  point  on 
the  elastic  axis.  Since  the  result  of  the  experimental 
tests  was  a  2  x  2  flexibility  influence  coefficient  matrix, 
the  same  result  is  desired  here,  to  allow  direct  comparison. 
The  2x2  matrix  requirement  dictates  that  we  use  the  equa¬ 
tions  of  motion  derived  from  the  two  term  deflection  equa¬ 
tion.  Observing  that  acceleration  (q)  for  static  deflec¬ 
tions  is  zero,  Eqs.  2-25a  and  b  are  repeated  here  with  q 
terms  eliminated. 


^11C 

^1  ^  ~ ip  ^7  ^  +  ^2  ^ 


2D 


16 

T~ 


1  -  Qi 


( 2-30a) 


2D 


*1  [ 


16  T 
7~  I< 


1  +  q,  ( 


CD 
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4D 


122. 


J11  + 


66 


cJl 


h  1  *  °2 


(2-30b) 
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The  generalized  forces  (Q^)  obtained  from  Eq.  2-20  are 


Q1  =  W  IZ  ] 


Q2  =  Va(  )/C 


( 2-31a) 
( 2-31b) 


where  Ffc  is  a  concentrated  test  load  and  M  is  a  concentrated 
test  moment  such  that 

£  -  Z  - 

F  =  /  f2  p  dydx  and  M  -  /  f2  y  p  dydx 

T  0  2  T  0  Z 


The  test  forces  and  moments  are  applied  at  the  same  loca¬ 
tions  and  with  the  same  sign  conventions  as  F£  and  in 

Fig.  2b.  By  the  relations  of  Eqs .  2-9  and  2-10,  the  def¬ 
lections  w_  and  a  taken  at  x  =  %Z  are 
E  4 

w£  =  )  qx  (2-32) 

ct  -  [  4>a  (  )/c  ]  q2  (2-33) 


Assuming  unit  forces  and  moments,  instead  of  unit 
displacements,  dictates  that  one  will  obtain  flexibility 
inflence  coefficients  (c^)  by  solving  the  equations  of 
motion.  The  matrix  of  influence  coefficients  must  then  be 
inverted  to  obtain  the  stiffness  influence  coefficients 
(k„)  .  This  is  generally  much  easier  to  accomplish  experi¬ 
mentally  than  to  try  and  obtain  the  stiffness  influence 
coefficients  directly,  by  applying  unit  displacements  to 
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the  structure.  To  aid  in  conceptualizing  the  problem,  the 
equations  of  motion  in  terms  of  flexibility  influences  are 
presented. 


*  C11FT  +  C12MT 


( 2-34a) 


a  *  c21FT  +  C22MT 


(2-34b) 


The  flexibility  influence  coefficients  are  defined  as  follows. 


°11E  Vft 


M  =  0 
T 


C12H  WE/MT 


F  =  0 
T 


(2-35) 


C21=-  “"'Iv  0 


c22=  a/Mtp 


V  0 


To  obtain  the  four  flexibility  influence  coefficients, 
Eqs.  2-30a  and  b  must  be  solved  simultaneously  twice,  each 
time  with  a  different  combination  of  and  Q2*  The  first 
combination  (FT  =  1  and  =  0)  will  be  called  "Test  1", 
and  the  second  combination  (FT  =  0  and  MT  =  1)  will  be  called 
"Test  2".  The  values  for  q^  and  q2  obtained  from  Test  1, 
when  substituted  into  Eqs.  2-32  and  2-33,  yield  c^  and  C2 ^ 
as  defined  by  Eq.  2-35.  Similarly,  q^  and  q2  from  Test  2 
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yield  c^2  and  c22.  The  final  step,  if  desired,  is  to 
arrange  the  flexibility  influence  coefficients  in  a  2  x  2 
matrix  and  invert  it  to  obtain  the  stiffness  influence  co¬ 
efficients  . 


2 . 4  Aeroelastic  Divergence  Problem 


Aeroelastic  divergence  is  a  static  deflection  problem 
such  that  the  aerodynamic  forces  and  moments  applied  to  the 
cantilevered  plate  are  a  function  of  the  torsional  deflec¬ 
tion  of  the  plate.  Divergence  is  defined  as  the  point  at 
which  the  restoring  forces  generated  by  the  structure  can 
no  longer  counteract  the  aerodynamic  forces.  At  that  point, 
by  the  simple  linear  theory  used  in  this  study,  the  deflec¬ 
tions  increase  without  bound  resulting  in  structural  failure 
Actually,  as  discussed  in  Chapter  III,  both  structural  and 
aerodyanmic  nonlinearities  come  into  play  at  large  deflec¬ 
tions  and  limit  the  maximum  deflections  to  finite  values. 
However,  structural  failure  may  still  occur.  These  non- 
linearities,  while  important,  will  not  be  considered  in  this 
development. 

The  static  forms  (q  =  0)  of  Eqs .  2-25a  and  b  are  re¬ 
peated  here. 


+  ^7 


Q,  ( 2-36a) 
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66 


cl 


^11° 

(  T  +  -  I  )  ] 

48D66l! 


*  Q. 


(2-36b) 


The  generalized  forces  (Q^)  obtained  from  Eq.  2-20  are 


/  4>1Le  dx 


( 2-37a) 


l  4> 


'  IT  “e 


dx 


{2—3  7b ) 


where 


L 


E 


“e 


pz  dy 
Y  Pz  dy 


L_  is  the  lift  per  unit  length  acting  at  the  elastic  axis 
of  the  plate,  in  the  same  direction  as  FE  (Fig.  2b) . 
is  the  aerodynamic  moment  per  unit  length,  again  acting 
at  the  elastic  axis. 

For  this  simple  static  analysis,  two-dimensional 
aerodynamic  strip  theory  was  deemed  to  provide  sufficient 
accuracy.  The  lift  per  unit  length  and  moment  per  unit 
length  are 

Le  =  qcaQ(a0  +  ae)  (2-38a) 


Mj.  =  qceaQ(a0  +  a&) 


( 2-38b) 
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where  q  is  the  dynamic  pressure,  e  is  the  distance  between 
the  elastic  axis  and  the  aerodynamic  center  (quarterchord) 
of  the  wing,  aQ  is  the  two  dimensional  lift  curve  slope, 
otg  is  the  rigid  plate  angle  of  attack  and  ag  is  the  flexible 
twist  of  the  plate,  such  that  an  +  a.  =  a  (Fig.  2b) . 

The  lift  curve  slope  of  a  flat  plate  of  infinite 
length  (span)  is  the  well  known  value  2 it  .  This  two  dimen¬ 
sional  lift  curve  slope  can  be  empirically  corrected  to 
account  for  the  affects  of  finite  span  by  applying  the 
correction  suggested  in  Ref.  13. 


3C, 


3a 


.  a  [ 

0l  AR  +  2 


(2-39) 


where  3C^/3a  is  the  corrected  lift  curve  slope  and  AR  is 
the  aspect  ratio  (2£/c)  of  an  equivalent  wing,  as  defined 
for  an  aircraft  wing  composed  of  two  cantilevered  plates. 

Equations  2-38a  and  b  must  be  transformed  to  generalized 
coordinates  before  they  are  incorporated  in  the  equations 
of  motion.  As  before,  Eq.  2-10  is  used  to  transform  aQ 
only.  The  generalized  forces  finally  take  the  form: 


Qx  =  qc2. 
Q2  =  qei 


5 

3a 

3a 


[  a0ri  +  c  I3  1 


[  aQI2  + 


c  X5  5 


(2-40a) 


( 2-40b) 


where  the  nondimensional  integral  expressions  are 


1 
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<f1  dx 

♦c.  dx 
*i+c  dx 
♦a  dx 


Substituting  Eqs .  2-40a  and  b  into  Eqs .  2-36a  and  b  yields 
the  complete  equations  of  motion.  These  equations  are 
presented  here  in  matrix  form  to  make  the  solution  technique 
more  apparent. 
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3a 
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3a 


)  ar 


(2-41) 


Classical  divergence  is  normally  evaluated  at  zero 


initial  angle  of  attack  (a^  =  0) ,  which  is  also  assumed 
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here.  Therefore,  a  unique  solution  to  Eq.  2-41,  with 
aQ  =  0,  can  be  found  by  setting  the  determinant  of  the  coef¬ 
ficient  matrix  equal  to  zero  and  solving  for  the  dynamic 
pressure  (q) .  The  dynamic  pressure  which  makes  the  deter¬ 
minant  go  to  zero  is  the  divergence  dynamic  pressure  (qD) , 
which  is  related  to  the  divergence  velocity  (UD)  by 


After  evaluating  the  determinant,  we  obtain  the  following 
expression  for  qD. 
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(2-43) 


When  the  integrals  are  evaluated  using  the  beam  mode 
shapes  in  Appendix  B,  the  contribution  of 
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will  be  small.  The  contribution  of 


DUD66 


to  the  value  of  q_,  while  not  negligibly  small,  is  much  less 
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than  that  of 
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The  denominator  of  Eq.  2-44  will  be  zero  when 
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’ll 
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2tl3I6 


=  0 . 3  —  (for  examples  in 

this  study) 


(2-44) 


When  equals  this  critical  value,  qQ  becomes  infinite 

and  the  wing  will  not  diverge.  Above  the  critical  value, 
qQ  is  negative  and  no  real  divergence  occurs  because  UD  will 
be  an  imaginary  number.  Below  the  critical  value,  including 
all  negative  values  for  D-^/D.,  3/  there  will  be  a  finite, 
real  UD. 

The  physical  mechanisms  which  cause  this  to  be  so  are 
readily  apparent  in  Fig.  3a  and  b.  Figure  3a  is  a  plate 
with  a  substantial  positive  D^.  In  this  case,  aerodynamic 
forces  arising  from  plate  displacements  tend  to  return  it 
to  equilibrium,  while  just  the  opposite  occurs  with  the  plate 
having  negative  Dlg  (Fig.  3b) . 

For  the  case  where  there  is  a  known,  nonzero,  initial 
angle  of  attack,  and  a  known  dynamic  pressure,  the  static 
deflections  of  the  structure  (a  and  wE)  can  be  obtained  by 
solving  the  equations  of  Eq.  2-41  simultaneously.  For  this 
study,  only  tip  deflections  are  evaluated,  because  they 
were  obtained  experimentally. 
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Eq.  2-41  can  be  rewritten  as 


kllql  +  k12q2  =  ala( 


(2-45a) 


k21ql  +  k22q2  "  a2a( 


( 2-45b) 


where 


Dllc 


2D...  ac. 

16  „  l  .. 

~V~  T6  "  qi  3a  1 3 


f  4D66  (i  ,  DHC  T 

[  cl  {18  JW77Tr  Iil 

66 


,  3C. 

gel  _ l_  , 

c  3a  15 


qc*  aT1  xi 


qei  3a  I2 


By  solving  Eqs .  2-45a  and  b  simultaneously,  one  obtains  the 
following  expressions  for  q^  and  q^. 


a2a0  ~  k22q2 


( 2-46a) 


[  aLk21  ~  *2kll  1 
[k12k21  “  kllk22J 


( 2-46b) 


The  generalized  deflections  are  then  transformed  to 
actual  deflections  (a  and  w_)  using  Eqs.  2-9  and  2-10, 
where  4> ^  and  <$>a  are  evaluated  at  x  =  l.  As  before,  a 
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(Fig.  2b)  equals  the  sum  of  cXq  and  .  The  total  angle  of 
attack  and  bending  deflection  at  the  tip  are 


<J>  U=£) 

aTIP  “  t  1  +  ( 


alk21  ~  a2kll 
k12k21  “  kllk22 
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(2-47a) 
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2.5  Free  Vibration  Problem 


The  free  vibration  problem  is  formulated  by  setting 
the  generalized  forces  (Q^)  equal  to  zero  in  either  Eqs . 
2-25a  and  b,  or  Eqs.  2-29a,  b  and  c.  To  make  the  equations 
of  motion  into  algebraic  equations  instead  of  differential 
equations,  one  assumes  that  the  motions  of  the  plate  will 
be  harmonic.  For  harmonic  (sinusoidal)  motion,  the  gener¬ 
alized  displacements  can  be  expressed  as 


q  =  qe 


iwt 


«  -  2,iwt 
q  =  -u  qe 


(2-48) 


where  co  is  the  frequency.  These  expressions  are  substituted 
into  the  differential  equations  of  motion  (Eqs.  2-25  and 
2-29)  to  obtain  the  sinusoidal  equations  of  motion.  The 
sinusoidal  equations  of  motion  are  presented  here,  in  matrix 


form,  for  both  the  two  term  and  three  term  deflection 
equations . 
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vibration  frequencies  (wn)  are  obtained  by  taking  the  square 
root  of  each  eigenvalue.  From  Eq.  2-49,  one  obtains  two 
natural  frequencies  (first  bending  and  first  torsion),  and 
from  Eq.  2-50,  one  obtains  three  natural  frequencies  (first 
bending,  second  bending  and  first  torsion) . 

To  simplify  calculations,  the  warping  stiffness  term 


48D66*2 


has  been  eliminated  because  its  contribution  is  small  for 
plates  with  moderate  to  high  aspect  ratios.  The  frequencies 
obtained  using  the  mode  shapes  in  Appendix  B  are  compared 
in  Chapter  IV  with  experimentally  obtained  frequencies. 

2 . 6  Aeroelastic  Flutter  Problem 

The  flutter  problem  is  formulated  using  an  indirect 
method  widely  known  as  the  U-g  method.  In  this  method, 
the  structural  damping  coefficient  (g) ,  introduced  into 
the  equations  of  motion,  is  plotted  versus  velocity  for 
each  vibration  mode.  Since  solutions  to  the  equations  of 
motion  represent  conditions  for  neutral  stability,  the 
value  of  g  obtained  in  this  manner  represents  the  amount 
of  damping  that  must  be  added  to  the  structure  to  attain 
neutral  stability  (flutter)  at  the  given  velocity.  Therefore, 


negative  values  of  structural  damping  indicate  that  the 
structure  is  stable.  Flutter  will  occur  when  the  artificial 
structural  damping  equals  the  actual  damping  of  the  structure 
To  simplify  calculations,  the  flutter  problem  will  be  formu¬ 
lated  using  only  the  equations  of  motion  from  the  two  term 
deflection  equation. 

The  equations  of  motion  {Eqs.  2-24a  and  b)  ,  as  derived 
in  Section  2.2,  again  assuming  sinusoidal  motion  and  neg¬ 
lecting  the  warping  stiffness  contribution,  are 
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and  the  generalized  forces  (Q.)  obtained  from  Eq .  2-20  are 
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(2-52b) 


Since  aerodynamic  forces  and  moments  obtained  by  assuming 
steady  flow  are  generally  inadequate  for  realistic  flutter 
calculations,  complete  unsteady  aerodynamic  expressions 
will  be  used.  Unsteady  aerodynamic  lift  (Lg)  and  moment 
(Mg)  expressions  derived  in  Ref.  13  have  been  suitably  trans 
formed  to  the  coordinate  system  in  Fig.  2.  The  reader  is 
cautioned  that  the  coordinate  system  used  throughout  this 
study  differs  from  the  standard  aeroelastic  coordinate 
system  used  in  Ref.  13.  The  difference  being  the  +wE  here 
is  -h  in  Ref.  13.  The  aerodyanmic  expressions  are 

Lg  =  7Tpb2[-Wg  +  Ua  -  baa  ]  +  2;rpUbC  (k)  [-w£ 

+  Ua  +  b ( 1/2  -  a)  a  ]  (2-53a) 

Mg  =  Trpb3  [~aWg  -  U(l/2  -  a)a  -  b(l/8  +  a2)a]  + 

rrpUb2  ( 1/2  +  a)  C(k)  [-Wg  +  Ua  +  b(l/2  -  a)  a  1 


(2-53b) 


43 


where  b 
a 

k 

C(k) 


semichord 

(distance  EA  is  aft  of  midchord) /b  (zero  in 
this  study) 

reduced  frequency  (wb/U) 

Theoaorsen  function 


Assuming  sinusoidal  motion 
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the  aerodynamic  expressions,  after  some  algebraic  manipu¬ 
lation,  become 
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[  M3  +  iM4  ]  =  (1/8  +  a2)  +  -(1-  2^-  C (k) 

+  i-  [(1/2  -  2a2)  C  (k)  -  (1/2  -  a)  ] 


Equations  2-54a  and  b  are  substituted  into  Eqs .  2-52a  and  b. 
One  observes  that,  since  the  plate  in  this  study  has  a  con¬ 
stant  chord,  variables  b,  [  +  il^  ),  [  +  iL^  ], 

[  +  iM2  ] ,  and  [  +  iM4  ]  are  not  functions  of  x  and 

can  be  brought  outside  the  integral.  Making  use  of  the  trans¬ 
formation  equations  (2-9  and  2-10) ,  and  expressing  the  integ¬ 
rals  in  nondimens ional  form,  we  arrive  at  the  final  expres¬ 
sions  for  the  generalized  forces. 
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Substituting  Eqs .  2-56a  and  b  into  Eq.  2-51  and  can¬ 
celing  elut  yields  a  new  form  of  the  equations  of  motion, 
which  are  written  here  in  contracted  matrix  form. 
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where  the  elements  of  the  K. .  and  A. .  matrices  are  given  as 
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Structural  damping  is  proportional  to,  and  opposes, 
the  deflection  of  the  structure.  The  structural  damping 
coefficient  (g)  is  introduced  into  Eq.  2-56  by  multiplying 
the  B  matrix  by  (1  +  ig) .  The  methodology  behind  this  step 
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is  treated  in  depth  in  Ref.  13.  An  assumption  here  is  that 
the  linear  damping  coefficient  (g^)  is  equal  to  the  rotary 
damping  coefficient  (g^) .  It  is  also  convenient  to  divide 
Eq.  2-56  by  -w2  and  combine  it  with  the  (1  +  ig)  term  to 
form  a  complex  eigenvalue  (Z) ,  defined  as 

Z  5  (  -+--i-g  )  (2-57) 

u>2 

Equation  2-53  now  becomes 

[  A  -  KZ  ]  q  =0  (2-58) 

where  the  elements  of  g  are  real  and  the  elements  of  4  are 
complex.  The  matrix  [  4  -  KZ  3  is  put  into  standard  eigen¬ 
value  form  by  premultiplying  by  K-1,  resulting  in 

[  K-1A  -  IZ  ]  2  =  0  (2-59) 

where  I  is  an  identity  matrix.  A  unique  solution  to  Eq. 

2-59  is  obtained  by  setting  the  determinant  of  the  coeffi¬ 
cient  matrix  equal  to  zero,  and  solving  for  Z.  For  a  2  x  2 
matrix,  the  two  Z- values  correspond  to  the  two  plate  vibra¬ 
tion  modes  (first  bending  and  first  torsion) . 

The  solution  procedure  is  to  pick  a  value  for  the  re¬ 
duced  frequency  (k)  and  solve  Eq.  2-59  for  Z^.  The  frequency 
structural  damping  coefficient  and  velocity  are  obtained 
from  Z  by  the  following  relations. 
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u  =  — - —  g  =  — { U  *  bw/k  (2-60a,b,c) 

/EeTTF  Re{  Z } 

A  new  k  value  is  chosen  and  the  procedure  is  repeated. 

Finally,  a  U-g  diagram  is  created  by  plotting  the  structural 
damping  for  each  root  (Z)  versus  velocity,  or  a  nondimensional 
flutter  speed  (U/bco^)  ,  were  uia  is  the  torsional  natural  vib¬ 
ration  frequency.  One  can  get  a  reasonable  approximation  for 

to  by  picking  a  very  large  value  for  k  and  solving  Eqs .  2-59 
ct 

and  2-60a  for  the  torsion  root.  This  is  essentially  the 
structural  vibration  frequency  when  U  equals  zero.  U-g 
diagrams  for  all  the  laminated  plates  considered  in  this 
study  appear  in  Appendix  B.  The  flutter  velocity  (Up)  was 
conservatively  chosen  to  be  the  point  where  the  damping 
coefficient  from  either  root  first  crosses  the  U-axis. 

Solving  the  flutter  problem  is  further  complicated  by 
the  presence  of  the  Theodorsen  function  (C(k)).  The  Theo- 
dorsen  function  is  a  nonrational  functional  which  can  be 
expressed  exactly  in  terms  of  modified  Bessel  functions  of 
the  second  kind.  Since  Eq.  2-59  was  to  be  implemented  on 
a  digital  computer  and  solved  many  times  for  different 
values  of  k,  an  approximation  for  C(k)  was  desired.  By 
taking  the  Laplace  Transform  of  R.T.  Jones'  exponential 
approximation,  one  obtains  a  rational  expression  for  C(k). 

This  expression  is16 
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0.5s2  +  0.2808s  +  0.01365 

C  { s )  =  - = -  (2-61) 

+  0.3455s  +  0.01365 

where  s  is  ik. 

The  flutter  problem  is  now  completely  defined.  For 
actual  computation,  taking  the  determinant  of  the  2x2 
coefficient  matrix  is  straightforward,  the  only  complica¬ 
tions  being  that  (1)  the  K-matrix  must  first  be  inverted 
and  (2)  that  the  coefficients  of  the  A-matrix  are  all  com¬ 
plex.  Because  of  the  simplicity  of  the  problem,  the  temp¬ 
tation  is  to  solve  it  manually.  However,  be  advised  that  a 
typical  U-g  plot  requires  on  the  order  of  20  data  points. 
Each  data  point  requires  that  Eq.  2-59  be  solved  for  a  new 
k  value. 
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CHAPTER  III 
EXPERIMENTS 


3.1  Test  Specimen  Selection 

Prior  to  making  any  test  specimens,  criteria  which  de¬ 
fined  desirable,  and  sometimes  essential,  characteristics 
of  the  test  specimens  were  established.  These  criteria  in¬ 
cluded 

(1)  Test  specimens  must  exhibit  a  wide  range  of  bending- 
torsion  coupling  stiffness. 

(2)  Test  specimens  must  be  rectangular,  constant  thickness, 
flat  plates  with  zero  sweep.  The  requirement  for  a 
flat  plate  dictated  using  a  symmetric  laminate,  since 
an  unsymmetric  laminate  would  warp  during  cure. 

(3)  Test  specimens  would  be  made  from  Hercules  AS/3501-6 
graphite/epoxy  since  it  was  available  at  M.I.T. 

(4)  Test  specimens  would  have  an  aspect  ratio  of  four, 
since  this  would  be  representative  of  an  aircraft  wing 
with  an  aspect  ratio  of  eight. 

(5)  Test  specimens  should  exhibit  flutter  and  divergence 
within  the  0-30  m/sec  speed  range  of  the  acoustic 
wind  tunnel  at  M.I.T. 

(6)  The  test  specimens  should  be  small  enought  to  be  made 
using  standard  TELAC  (Technology  Laboratory  for  Ad¬ 
vanced  Composites)  ply  cutting  templates  and  curing 
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plates.  Additionally,  a  small  size  was  desirable  be¬ 
cause  the  effect  of  slight  material  warping,  often 
encountered  during  curing,  would  be  minimized. 

(7)  The  laminate  should  be  able  to  withstand  repeated 
large  static  and  oscillatory  loads. 

Criteria  five  and  six  directly  oppose  one  another  when 
trying  to  arrive  at  an  optimum  design.  Trying  to  satisfy 
both  of  them  dictated  a  very  thin  plate.  Also,  to  achieve 
substantial  bending- tors ion  coupling,  an  unbalanced  lamina¬ 
tion  scheme  was  chosen.  Unbalanced  laminates  have  an  unequal 
number  of  +0  and  -6  plies.  However,  a  midplane  symmetric 
ply  arrangement  was  still  required.  Carrying  the  unbalance 
to  the  limit  results  in  a  unidrectional  laminate  with  all 
plies  at  +0  or  all  plies  at  -9 .  To  improve  the  toughness 
(criterion  seven)  of  a  unidirectional,  off-axis  laminate,  the 
center  two  plies  would  have  a  zero-degree  ply  angle.  The 
total  thickness  was  chosen  to  be  six  plies.  Five  different 
laminates  were  selected  as  representative  of  a  wide  range 
of  bending-torsion  coupling  stiffnesses,  both  positive  and 
negative.  The  first  two  laminates,  [+30-/0]  and  [-30-/0]  , 

m  S  £  9 

were  chosen  because  a  unidirectional  laminate  with  0  =  30° 
exhibits  the  highest  bending-torsion  coupling  stiffness 
(D^g)^.  The  next  two  laminates,  t-+45  ]  s  ana  [-452/0]  g, 

were  chosen  because  a  unidirectional  laminate  with  8  =  45° 
exhibits  the  highest  torsion  stiffness  (D--) ,  but  still  has 
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a  large  D^g.  The  final  laminate,  [0^/90] s,  was  chosen  be¬ 
cause  a  unidirectional  laminate  with  0=0°  exhibits  the 
highest  bending  stiffness  (D^)  ar*d  2ero  Di6*  During  the 
test  program,  another  laminate,  [*45/0]  ,  was  added  because 
it  had  the  same  theoretical  and  Dgg  as  [  +  452/0]  ,  with 
less  D^g .  To  minimize  the  number  of  test  specimens  which 
had  to  be  constructed,  the  [+82/0] s  and  [-92/0]s  test  spe¬ 
cimens  were  actually  the  same  laminate,  simply  rotated 
180  degrees  about  the  x-axis  (Fig.  1) 

The  test  specimens  would  have  to  have  an  overall 
length  of  330  mm  (13  in)  and  a  chord  of  76  mm  (3  in) ,  be¬ 
cause  this  was  the  largest  size  that  could  be  constructed 
satisfying  criteria  four  and  six.  The  overall  length  in¬ 
cluded  a  25  mm  (1  in)  loading  tab,  making  the  effective 
cantilever  plate  length  305  mm  (12  in) .  Finally,  in  an 
effort  to  minimize  plate  stiffness,  an  airfoil  shaped 
fairing  would  not  be  used  on  the  plate  for  wind  tunnel 
tests.  This  had  the  added  advantage  of  making  the  plate 
stiffness  easier  to  calculate. 

3 . 2  Test  Specimen  Preparation 

The  test  specimens  were  constructed  from  Hercules 
AS/3501-6  graphite/epoxy  prepreg  tape  from  Lot  No.  1643. 

The  tape  was  305  mm  (12  in)  wide  and  had  a  nominal  thick¬ 
ness  of  0.134  mm  (0.00528  in).  Individual  plies  were  cut 
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to  the  proper  size  and  angular  orientation  using  aluminum 
templates,  and  assembled  into  305  mm  (12  in)  by  356  mm 
(14  in)  laminates.  The  laminates  and  curing  materials  (Ap¬ 
pendix  C)  were  arranged  on  an  aluminum  curing  plate  especi¬ 
ally  designed  by  students  in  TELAC  for  use  with  graphite/ 
epoxy  laminates  having  a  length  and  width  as  stated  above. 

The  laminates  were  cured  in  a  Baron  model  BAC-35  autoclave 
using  the  cure  cycle  listed  in  Appendix  C.  After  curing, 
the  laminates  were  post-cured  in  a  forced  air  circulation 
oven  at  350°  F  for  eight  hours.  After  post-curing,  a  rec¬ 
tangular  test  specimen  330  mm  (13  in)  long  and  76  mm  (3  in) 
wide  was  cut  from  each  of  the  four  laminates  using  a  diamond 
coated  cutting  wheel  mounted  on  an  automatic  feed,  milling 
machine. 

The  thickness  of  a  graphite/epoxy  plate  tended  to 
vary  over  its  surface.  Therefore,  the  plate  thickness  was 
measured  at  several  different  locations  and  averaged.  The 
same  procedure  was  followed  in  measuring  the  length  and 
width,  although  the  variation  was  much  less.  The  averaged 
measurements  for  each  test  specimen  appear  in  Appendix  D, 
along  with  the  nominal  values.  The  thickness  variation 
between  laminates  (all  were  six  ply)  was  0.033  mm  (0.0013  in) 
with  the  average  being  0.807  mm  (0.0318  in).  This  compared 
favorably  to  the  nominal  thickness  of  0.804  mm  (0.0317  in) 
for  a  six  ply  laminate.  The  laminates  were  also  weighed 
on  a  triple  beam  balance,  from  which  the  material  density 
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(p)for  each  test  specimen  was  calculated.  The  average  den- 

sity  was  1.53  x  10  kg/m  (0.0554  lb/in °)  with  a  maximum 

variation  of  0.06  x  103  kg/m3  (0.0022  lb/in3) .  This  again 

3  3 

compared  favorably  with  a  nominal  density  of  1.52  x  10  kg/m 
(0.055  lb/in^)  for  graphite  epoxy. 

Loading  tabs  83  mm  (3.25  in)  x  25.5  mm  (1  in)  were 
machined  from  2.4  mm  (0.094  in)  aluminum  plate  and  bonded  to 
the  base  of  each  test  specimen  with  epoxy,  cured  at  room 
temperature.  The  loading  tabs  were  intended  to  aid  in  align¬ 
ing  the  test  specimen  in  the  clamping  fixture  and  to  prevent 
damage  to  the  plate  surface  fibers . 

To  get  an  indication  of  the  lateral  deflections  the 
plate  would  experience  during  the  wind  tunnel  tests,  strain 
gauges  were  attached  to  the  base  of  each  test  specimen  at 
the  midchord,  as  shown  in  Fig.  4.  Two  Micro  Measurements 
EA-Q3-1873B-120  strain  gauges,  from  Lot.  No.  R-A21AB02  with 
a  gauge  factor  of  2.085,  were  attached  to  each  test  specimen 
(one  on  each  side)  to  measure  bending  strain.  Two  BLH-SR4 ; 
FAED-25B-12-59 ;  Serial  No.  5-AE-SC  strain  gauge  rosettes 
from  Lot  No.  A-315  with  a  gauge  factor  of  2.01  were  attached 
to  each  test  specimen  (one  on  each  side)  to  measure  torsion 
strain.  The  two  bending  gauges  were  wired  together  as  a 
two-arm  bridge  circuit  with  three  external  lead  wires  approxi¬ 
mately  305  mm  (12  in)  long.  The  two  torsion  gauges  were 
wired  together  as  a  four-arm  bridge  circuit  with  four  external 
lead  wires.  Wiring  the  strain  gauges  on  either  side  of  the 
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Figure  4.  Sample  Test  Specimen 
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plate  together,  in  this  manner,  doubled  the  signal  output 
for  each  channel  (bending  and  torsion)  and  provided  auto¬ 
matic  temperature  compensation.  This  temperature  compen¬ 
sation  was  very  important  during  the  wind  tunnel  tests,  where 
air  flowing  over  an  uncompensated  gauge  would  cause  a  zero 
shift  as  the  tunnel  velocity  was  varied.  The  final  step, 
after  all  solder  connections  had  been  made,  was  to  coat  the 
gauges  with  Micro  Measurements  M-Coat  A,  an  air-drying  poly¬ 
urethane,  for  protection. 

3 . 3  Test  Apparatus  and  Procedure 
(a)  Static  Deflection  Tests 

The  static  deflection  test  setup  is  shown  in  Fig.  5. 

It  consisted  of  a  330  mm  (13  in)  x  508  mm  (20  in)  plywood 
base  with  six  vertical  steel  rods  approximately  760  mm 
(30  in)  long.  The  test  specimen  was  clamped  in  a  vise 
machined  from  a  25.4  mm  (1  in)  x  152  mm  (6  in)  x  229  mm 
(9  in)  aluminum  block,  which  was  bolted  to  the  base  of  the 
test  fixture.  Two  removable,  low  friction,  pulleys  were 
attached  to  the  vertical  rods  such  that  a  force  or  moment 
could  be  applied  to  the  test  specimen  at  any  location  along 
its  length.  Rulers,  graduated  in  32nd 's  of  an  inch,  were 
also  attached  to  the  vertical  rods  to  facilitate  measuring 
the  test  specimen's  edge  deflections.  A  deflection  indicator 
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was  constructed  out  of  balsa  wood.  It  had  needle  pointers 
on  the  ends  and  cotton  threads  attached  to  its  midpoint  and 
ends.  The  threads,  when  routed  over  the  pulleys  and  attached 
to  weights,  transferred  a  force  or  moment  to  the  test  speci¬ 
men  (Fig.  5)  . 

The  deflection  indicator  was  aligned  with  the  lines 
scribed  on  the  test  specimen  at  the  3/4  length  point,  and 
the  test  specimen  was  clamped  in  the  vise.  The  pulleys  were 
clamped  to  the  middle  two  vertical  rods  (Fig.  5a)  at  the 
proper  height,  and  threads  from  the  center  of  the  deflec¬ 
tion  indicator  were  routed  over  the  pulleys.  The  rulers 
were  adjusted  to  the  proper  height  and  zeroed  with  respect 
to  the  deflection  indicator  pointers.  Weights  of  10,  20, 

30,  and  40  grams  were  successively  attached  to  the  threads, 
first  to  give  positive  deflections,  then  to  give  negative 
deflections.  As  each  weight  was  attached,  the  readings 
from  both  pointers  were  recorded,  along  with  the  weight. 

Next,  the  pulleys  were  moved  to  the  corner  rods  (Fig. 

5b)  and  the  end  threads  were  routed  over  them,  so  as  to  pro¬ 
vide  a  positive  moment  when  the  weights  were  attached.  Weights 
of  10,  20,  33.9,  and  43.9  grams  were  successively  attached  to 
each  thread  of  the  couple,  and  readings  from  both  pointers 
were  again  recorded  along  with  the  weights.  The  pulleys  were 
then  switched  to  the  corner  rods  on  the  opposite  diagonal, 
and  the  procedure  was  repeated  for  a  negative  moment. 
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For  each  data  point,  the  lateral  deflection  of  the 
elastic  axis  (w£)  and  the  rotation  of  the  test  specimen 
about  the  elastic  axis  (a)  were  calculated  using  the  data 
reduction  formulas  in  Appendix  E.  The  lateral  deflections 
obtained  from  the  load  test  (Test  1)  were  plotted  versus 
load,  with  the  slope  being  the  bending  flexibility  influence 
coefficient  (c^)  .  The  angular  deflections  (a)  from  Test  1 
were  plotted  versus  load,  yielding  the  bending-torsion 
coupling  flexibility  influence  coefficient  (c2^) •  The  an¬ 
gular  deflections  obtained  from  the  moment  test  (Test  2) 
were  plotted  versus  moment,  yielding  the  torsion  flexibility 
influence  coefficient  (c22) •  The  lateral  deflections  ob¬ 
tained  from  Test  2  were  plotted  versus  moment  to  obtain 
the  other  bending- tors  ion  coupling  flexibility  influence 
coefficient  (c12)  •  The  flexibility  influence  coefficient 
plots  for  the  four  test  specimens  appear  in  Appendix  E. 

The  flexibility  influence  coefficients  for  a  [+02/O]s  anc* 
a  [-02/O]g  laminate  will  be  the  same,  except  for  the  sign 
on  c12  and  c2^,  since  they  were  physically  the  same  test 
specimen.  Finally,  the  flexibility  influence  coefficients 
for  each  test  specimen  were  arranged  in  a  2  x  2  matrix  (C) 
and  inverted  to  obtain  the  stiffness  influence  coefficient 
matrix  (K) .  The  results  of  the  static  deflection  tests  are 
discussed  in  Section  3.4. 
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(b)  Free  Vibration  Tests 

The  free  vibration  test  setup,  shown  in  Fig.  6,  used 
the  same  vise  as  the  static  deflection  tests.  The  vice  was 
suspended  by  four  spring  steel  strips,  which  allowed  it  to 
translate  along  the  z-axis  (Fig.  2) ,  but  restricted  motion 
in  all  other  directions.  The  vise  was  rigidly  attached  to 
a  horizontally  mounted  Ling  model  420  shaker  using  a  cylin¬ 
drical  aluminum  adapter.  The  shaker  had  a  peak  force  of 
445  N  (100  lbs)  and  a  frequency  range  of  approximately  5  to 
3000  Hz.  Since  the  shaker  was  driven  by  an  audio  amplifier, 
there  was  a  certain  amount  of  distortion  to  the  sine  wave 
output  signal  below  20  Hz.  For  some  vibration  tests,  an 
Endevco  model  7701-50  "Isoshear"  accelerometer  was  mounted 
to  the  vise  by  a  threaded  mounting  stud,  and  an  Endevco 
model  222B  "Micro-miniature"  accelerometer  was  mounted  to 
the  test  specimen,  within  25  mm  of  the  base,  using  Eastman 
910  adhesive.  The  outputs  of  these  accelerometers,  after 
passing  through  special  amplifiers,  were  connected  to  a 
Tektronix  type  502  dual  beam  oscilloscope  so  that  both  out¬ 
puts  could  be  displayed  simultaneously.  Finally,  a  digital 
signal  counter  was  attached  directly  to  the  signal  generator 
to  provide  an  accurate  frequency  readout. 

The  test  specimen  was  aligned  and  clamped  in  the  vise. 

A  white  paper  screen  was  suspended  behind  the  test  specimen 
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Figure  6. 


Free  Vibration  Test  Apparatus 
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to  aid  in  observing  the  vibration  modes.  The  test  began 
with  the  signal  generator  set  at,  or  below,  3  Hz.  The  fre¬ 
quency  was  slowly  increased  until  the  first  bending  mode  was 
excited.  The  frequency  was  very  finely  adjusted  until  the 
vibration  amplitude  peaked,  as  observed  visually,  and  the 
frequency  was  recorded.  This  procedure  was  repeated  for 
the  second  bending  and  first  torsion  vibration  modes.  Since 
flutter  is  normally  associated  with  the  lower  (bending  and 
torsion)  vibration  modes ,  only  the  first  bending,  second 
bending  and  first  torsion  frequencies  were  obtained  for  each 
specimen. 

In  an  effort  to  improve  the  accuracy  and  repeatability 
of  the  frequencies,  the  dual  beam  oscilloscope  was  used  to 
observe  the  accelerations  of  the  vise  and  the  test  specimen. 
The  resonance  frequency  was  defined  as  the  frequency  where 
a  90°  phase  shift  occurred  between  the  two  accelerometer 
traces.  If  a  complete  90°  phase  shift  did  not  occur,  then 
the  resonance  frequency  was  taken  to  be  the  frequency  where 
the  maximum  phase  shift  occurred.  In  an  attempt  to  obtain 
a  clear  phase  shift  at  the  first  torsion  natural  mode,  the 
specimen-mounted  accelerometer  was  moved  to  different  chord- 
wise  locations,  all  the  way  to  the  edge  of  the  test  specimen 
The  results  of  the  free  vibration  tests  are  discussed  in 


Section  3.4. 
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(c)  Wind  Tunnel  Tests 

All  wind  tunnel  tests  were  performed  in  the  M.I.T., 
Department  of  Aeronautics  and  Astronautics,  acoustic  wind 
tunnel.  The  acoustic  wind  tunnel  is  a  continuous  flow  tun¬ 
nel  with  a  1.5  m  (5  ft)  x  2.3  m  (7.5  ft)  free  jet  test  sec¬ 
tion  2.3  m  (7.5  ft)  long,  located  inside  a  large  anechoic 
chamber.  The  tunnel  was  powered  by  a  100  HP  motor  giving 
it  a  continuously  variable  velocity  range  of  0  to  32  m/sec 
(0  to  105  ft/sec) .  The  tunnel  control  panel  was  located 
outside  the  anechoic  chamber.  The  velocity  was  controlled 
by  two  levers  (coarse  and  fine  speed  control) .  The  coarse 
lever  controlled  the  motor  field  current  and  was  variable 
in  fixed  step  increments  only.  The  fine  lever  controlled 
the  motor  shunt  current  and  was  continuously  variable. 

The  total  current  being  drawn  by  the  motor,  along  with 
voltage,  were  displayed  at  the  control  panel  on  analog 
gauges.  There  was  a  400  amp  current  limitation  on  the  motor, 
which  was  protected  by  circuit  breakers.  The  tunnel  velo¬ 
city  was  read  from  an  alcohol  manometer,  calibrated  in 
inches  of  alcohol,  and  located  at  the  control  panel.  The 
manometer  was  connected  to  a  pitot  tube  located  slightly 
forward  of  the  test  section.  A  hydraulic  motor,  recessed 
into  the  floor  of  the  test  section,  provided  a  means  of 
mounting  the  test  apparatus. 
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The  test  setup,  shown  in  Fig.  7,  consisted  of  a  turn¬ 
table  machined  from  aluminum  mounted  on  a  635  mm  (25  in) 
tall  pedestal  made  of  152  mm  (6  in)  steel  pipe.  The  pedes¬ 
tal  was,  in  turn,  mounted  to  the  floor  of  the  wind  tunnel 
test  section.  The  vise  from  the  static  deflection  tests 
was  fastened  to  the  free  rotating  portion  of  the  turntable. 

A  wooden  cover  disk  508  mm  (20  in)  in  diameter,  shown  re¬ 
moved  in  Fig.  7a  and  attached  in  Fig.  7b,  was  used  to  pro¬ 
vide  smooth  airflow  past  the  test  specimen.  The  disk,  which 
had  the  angle  of  attack  marked  on  its  edge  in  2°  increments 
up  to  18°,  rotated  with  the  test  specimen.  A  pointer  at¬ 
tached  to  the  fixed  base  provided  a  consistent  means  of 
reading  the  angle  of  attack.  An  aluminum  rod,  attached  to 
to  the  disk  and  extending  outside  the  test  section,  allowed 
the  angle  of  attack  to  be  varied  while  the  test  was  in 
progress. 

A  terminal  strip  attached  to  the  underside  of  the  disk 
provided  a  convenient  means  of  connecting  the  test  specimen 
strain  gauge  lead  wires  to  the  instrumentation  wiring  harness. 
The  instrumentation  system  was  a  Gould  2400  series  four  chan¬ 
nel  strip  chart  recorder.  The  recorder  had  a  four  arm  D.C. 
bridge  preamplifier  installed  in  channel  1  and  a  two  arm 
D.C.  bridge  preamplifier  installed  in  channel  2.  The  recorder 
also  had  an  internal  time  base  generator  which  marked  the 
strip  chart  at  one  second  intervals.  The  preamplif ier  had 
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the  capability  of  both  attenuating  and  amplifying  the  input 
signal.  The  strip  chart  speed  could  be  set  at  5,  10,  25,  50, 
100,  and  200  mm/sec  plus  each  speed  could  be  divided  by  100. 
The  chart  speed  was  accurate  to  *  0.5%.  To  provide  a  record 
of  the  test  specimen  deflections,  the  torsion  strain  gauge 
was  attached  to  channel  1  and  the  bending  strain  gauge  was 
attached  to  channel  2. 

The  test  specimen  was  aligned  and  clamped  in  the  vise. 
The  cover  disk  was  installed,  and  all  gaps  were  covered  with 
tape.  The  strain  gauge  lead  wires  were  attached  to  the 
terminal  strip,  and  the  strip  chart  recorder  was  adjusted 
according  to  the  manufacturer's  instructions.  The  test 
number  and  test  specimen,  along  with  test  section  tempera¬ 
tures  and  atmospheric  pressure,  were  recorded  on  a  data 
sheet.  The  first  three  natural  vibration  frequencies  were 
recorded  on  the  strip  chart  by  giving  the  test  specimen  a 
bending  deflection  and  letting  it  oscillate,  then  giving 
it  a  torsion  deflection  and  letting  it  oscillate.  Static 
calibrations  were  performed  by  giving  the  test  specimen  a 
series  of  known  pure  bending  and  torsion  tip  deflections 
and  appropriately  annotating  the  strip  chart  deflections. 

The  angle  of  attack  was  adjusted  to  zero,  and  the  wind 
tunnel  was  started  and  set  at  a  low  velocity  (<  8  m/sec; 

25  ft/sec) .  After  the  velocity  stabilized,  a  reading  was 
taken  from  the  alcohol  manometer  andannotated  on  the  strip 
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chart,  in  inches  of  alchol.  The  angle  of  attach  was  then 
swept  from  0°  to  12°  (sometimes  to  18°  if  deemed  appropri¬ 
ate)  in  2°  increments,  while  the  velocity  was  held  constant. 
Each  angle  of  attack  was  annotated  on  the  strip  chart.  The 
angle  of  attack  was  returned  to  zero  and  the  tunnel  speed 
increased  slightly.  The  new  manometer  reading  was  annotated 
on  the  strip  chart,  and  the  angle  of  attack  sweep  was  re¬ 
peated.  The  strip  chart  recorder  was  run  at  1  mm/sec  for 
most  of  the  test.  When  oscillation  was  detected,  a  two 
second  expanded  record  (100  or  200  mm/sec)  was  made.  For 
some  test  specimens,  at  higher  angles  of  attack,  the  oscil¬ 
lations  would  begin  very  abruptly.  To  get  better  data  on 
the  angle  of  attack,  where  flutter  actually  started,  the 
angle  of  attack  was  increased  in  1°  increments  in  this 
region.  This  test  procedure  was  repeated  until  either  the 
the  maximum  tunnel  speed  was  reached,  or  the  test  specimen 
deflections  became  excessive.  At  higher  tunnel  speeds,  the 
maximum  angle  of  attack  tested  was  reduced  to  prevent  damage 
to  the  test  specimen.  For  the  specimens  which  exhibited 
flutter  at  zero  angle  of  attack,  the  normal  procedure  was 
to  obtain  two  or  three  more  data  runs  at  velocities  slightly 
higher  than  the  velocity  where  the  test  specimen  fluttered 
at  zero  angle  of  attack. 

Divergence  testing  often  results  in  the  destruction  of 
a  test  specimen.  Since  this  was  not  a  desirable  result. 
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wind  tunnel  speeds  were  increased  in  very  small  increments 
as  the  divergence  velocity  was  approached.  The  deflections 
of  the  test  specimen  were  closely  monitored,  and  when  it 
would  no  longer  return  to  a  zero  deflection  state,  the 
tunnel  velocity  was  not  increased  further.  The  extreme 
flexibility  of  the  test  specimens  also  increased  the  allow¬ 
able  margin  for  error. 

Finally,  after  all  the  desired  data  was  taken,  the 
tunnel  was  shut  down  and  the  temperature  of  the  test  sec¬ 
tion  was  recorded  on  the  data  sheet.  Since  the  wind  tunnel 
was  a  continuous  flow  variety,  temperature  increases  of  up 
to  10°  F  were  encountered  during  test  runs.  The  test  speci¬ 
men  was  removed  from  the  vise  and  another  test  specimen 
installed.  When  going  from  a  [+0 2/°^  specimen  to  a 
[  —  0 2/0  3  s  specimen,  the  turntable  was  simply  rotated  180 
degrees,  while  the  specimen  remained  undisturbed.  The  test 
procedure  was  then  repeated  for  the  new  test  specimen. 

After  the  wind  tunnel  tests  were  complete,  the  mano¬ 
meter  readings  were  converted  to  indicated  velocity,  and 
then  to  true  velocity  using  the  conversions  in  Appendix  E. 
The  temperature  used  for  the  true  velocity  correction  was 
the  average  of  the  start  and  finish  temperature  for  each 
run.  The  oscillation  frequencies  were  obtained  manually 
and  annotated  on  the  strip  chart,  along  with  the  indicated 
velocity.  The  static  calibrations  were  used  to  obtain  an 
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approximate  relationship  between  test  specimen  deflection 
and  strain  gauge  output  using  the  procedure  outlined  in 
Appendix  E.  This  relationship  was  used,  on  selected  tests, 
to  calculate  the  test  static  deflections  (not  oscillation 
deflections) .  The  peak  to  peak  oscillation  amplitude  of 
the  bending  and  torsion  strain  gauges,  on  selected  tests, 
was  also  annotated  on  the  strip  chart.  Finally,  the  flutter 
boundary  velocities,  along  with  the  dimensionless  reduced 
flutter  speeds  (U/b<jJa)  were  plotted  versus  angle  of  attack 
for  each  test  specimen.  On  selected  tests,  the  test  speci¬ 
men  deflection,  for  a  given  angle  of  attack,  was  plotted 
versus  velocity.  Also,  flutter  amplitude,  for  a  given 
angle  of  attack,  was  plotted  versus  velocity  for  selected 
tests.  These  test  results  are  discussed  in  Section  3.4. 

3 . 4  last  Results 

(a)  Static  Deflection  Tests 

The  flexibility  influence  coefficients  (c^..)  which  are  a 
measure  of  the  structure's  deflection  for  a  given  load,  were 
obtained  by  determining  the  slope  of  the  lines  plotted  in 
Appendix  E.  All  the  plots  were  linear,  with  the  exception 
of  c22  f°r  the  [02/90]s  test  specimen.  For  this  plot,  a 
linear  approximation  was  used  to  determine  the  slope.  The 
data  scatter  on  all  the  plots  was  less  than  10%.  The  flexi¬ 
bility  influence  coefficients  could  have  been  arranged  in 
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a  2  x  2  matrix  and  inverted  to  obtain  the  stiffness  influence 
coefficients  (k^j) »  which  are  a  measure  of  the  force  a  struc¬ 
ture  exerts  for  a  given  displacement.  However,  for  this 
section,  the  flexibility  influence  coefficients  will  be  exa¬ 
mined  directly.  Table  1  lists  the  experimentally  obtained 
inflence  coefficients  for  each  test  specimen. 

TABLE  1 

EXPERIMENTAL  FLEXIBILITY  INFLUENCE  COEFFICIENTS 

Test  C11  c12  C21  c22  Coupling 

Specimen  (ft/lb)  (ft/ft-lb)  (rad/lb)  (rad/fflb)  Factor  (<) 


[02/9°  ]s 

0.18 

-0.03 

-0.02 

2.3 

0.0014 

[±45/0]g 

0.60 

-0.16 

-0.13 

1.06 

0.033 

[+452/°1s 

0.98 

-0.76 

-0.75 

2.7 

0.22 

[-452/°]s 

0.98 

+  0.76 

+  0.75 

2.7 

0.22 

[+302/0]s 

0.64 

-0.80 

-0.82 

3.0 

0.34 

[-3°2/°]s 

0.64 

+0.80 

+0.82 

3.0 

0.34 

A  measure  of 

the  accuracy 

of  these 

static  deflection 

tests  was  that,  according  to  Maxwell's  Law  of  Reciprocal 
Deflections,  c^2  must  equal  c21.  All  of  the  test  specimens 
showed  excellent  agreement  in  this  area. 

The  coupling  factor  (tc)  ,  included  in  Table  1,  was  de¬ 
fined  as 

_  C12C21 
<  =  - 

C11C22 
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and  was  a  dimensionless  measure  of  the  amount  of  bending- 
torsion  coupling  a  test  specimen  possessed.  One  observes 
that  going  from  a  balanced  to  an  unbalanced  laminate  sig¬ 
nificantly  increased  the  coupling.  Also,  the  [+302/0]g  and 
[-302/0]s  test  specimens  had  the  highest  coupling  of  those 
tested. 

In  going  from  the  balanced  [*45/0]  test  specimen  to 

S 

the  unbalanced  [+452/0]g  test  specimen,  the  torsional  flexi¬ 
bility  ( c 2 2 )  increased  by  a  factor  of  2.5,  and  the  bending 
flexibility  (c^)  increased  by  a  factor  of  1.6.  The  reader 
should  also  note  that  the  [+3O2/0]g  test  specimen  was  more 
flexible  in  torsion  than  the  [O2/90]g  test  specimen. 

(b)  Free  Vibration  Tests 

Excitation  of  the  torsion  vibration  mode  was  extremely 
difficult  on  the  [C>2/90]s  anci  [±45/0]gtest  specimens, 

which  had  very  little  bending- torsion  coupling.  Addition¬ 
ally,  the  torsion  mode  was  only  lightly  excited  on  the  other 
test  specimens.  The  accelerometer  bonded  to  the  test  speci¬ 
men  was  also  largely  unsuccessful  for  the  torsion  mode,  be¬ 
cause  the  phase  shift  was  very  slight.  Since  the  shaker 
was  driven  by  an  audio  amplifier,  the  output  signal  was 
distorted  below  20  Hz.  This  made  using  an  oscilloscope  to 
detect  phase  shift  in  the  first  bending  mode  difficult. 

The  phase  shift  technique  did,  however,  work  well  for  the 
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second  bending  mode . 

Frequencies  of  the  lowest  three  natural  vibration  modes, 
obtained  from  both  the  shaker  and  wind  tunnel  tests,  are 
tabulated  in  Appendix  E.  The  averaged  frequencies  for  each 
test  specimen  appear  in  Table  2.  Frequencies  from  both 
sources  were  all  within  3%  of  the  average.  The 
and  [ — 9 0 ] s  test  specimens,  being  the  same  physical  plate, 
had  the  same  natural  frequencies. 


TABLE  2 

EXPERIMENTAL  NATURAL  FREQUENCIES  (AVERAGE  VALUES) 


Test  1st  bending  2nd  bending  1st  torsion 

Specimen  freq  (Hz)  freq  (Hz)  freq  (Hz) 


[02/9°]s 

11.1 

69.0 

42.2 

[*45/0 ] s 

6.07 

38.5 

77.4 

[+45  /0]  ' 

) 

-VoJ 

>  4.83 

1 

30.1 

51.1 

[+30  /O], 

z  s  ( 

[-30-,/°]/ 

►  5.97 

1 

36.1 

58.4 

In  going  from  the  balanced  [*45/01  test  specimen  to 

5 

the  unbalanced  [+45-/0]  test  specimen,  the  torsion  fre- 

Cm  S 

quency,  which  is  largely  dependent  on  the  torsional  flexi¬ 
bility  of  the  specimen,  decreased  by  34%.  Also,  the  first 
bending  mode  frequency  decreased  by  20%.  These  results 
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were  consistent,  at  least  in  direction,  with  those  of  the 
static  deflection  tests. 

The  situation  becomes  slightly  more  confusing  when  the 
frequencies  of  the  other  test  specimens  are  compared  with 
their  flexibility  influence  coefficients.  The  first  and 
second  bending  mode  frequency  trends  correlate  well  with 
the  trends  of  the  bending  flexibility  influence  coefficients 
However,  the  trends  for  the  torsion  frequencies  and  the  tor¬ 
sion  flexibility  influence  coefficients  had  poor  correlation 
This  indicated  that  the  bending  torsion  coupling  must  be 
primarily  affecting  the  torsion  frequency.  As  an  aside, 
when  stiffness  influence  coefficients  were  compared  to 
natural  frequencies  in  this  manner,  one  was  led  to  conclude 
that  the  coupling  primarily  affected  the  bending  frequencies 

All  the  test  specimens,  except  the  [O2/90]s  specimen, 
exhibited  first  bending  as  the  lowest  frequency,  followed  by 
second  bending,  and  finally  first  torsion.  The  [02/ 9 0 ] s 
test  specimen  started  with  first  bending,  followed  by  first 
torsion  and  finally  second  bending. 

(c)  Wind  Tunnel  Tests 

The  primary  objective  of  the  wind  tunnel  tests  was  to 
determine  the  zero  initial  angle  of  attack  (<*q  =  0)  diver¬ 
gence  and  flutter  velocities  for  each  test  specimen.  Since 
the  test  velocities  were  low  and  the  test  specimens  quite 
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flexible,  the  onset  of  flutter  at  zero  aQ  tended  to  be  gra¬ 
dual.  For  consistency,  the  flutter  threshold  was  defined 
to  be  the  velocity  where  the  oscillation  amplitude  stabi¬ 
lized  and  the  waveform  became  periodic.  Determining  the 
threshold  for  stall  flutter  was  somewhat  easier,  because 
oscillation  tended  to  begin  abruptly  when  the  stall  angle 
of  attack  was  reached. 

The  type  of  flutter  (bending,  torsion  or  bending- 
torsion)  could  not  be  determined  completely  from  the  ampli¬ 
tude  of  the  bending  and  torsion  strip  chart  channels.  This 
was  because,  for  a  plate  with  bending  -  torsion  coupling,  the 
elastic  axis  was  not  located  at  the  midchord,  where  the 
strain  gauges  were  located.  So  a  pure  bending  deflection 
excited  both  the  bending  and  torsion  strain  gauges.  The 
problem  was  overcome  by  defining  bending  flutter  as  oscil¬ 
lations  with  a  frequency  very  close  to  one  of  the  bending 
natural  frequencies.  Torsion  flutter  was  said  to  occur 
when  the  frequency  was  close  to  the  torsion  natural  fre¬ 
quency,  and  bending-torsion  flutter  occured  at  a  frequency 
between  the  other  two. 

On  test  9,  with  the  1-452/0] g  test  specimen,  a  torsion 
strain  gauge  failed  due  to  excessive  deflection  about  half 
way  through  the  test.  The  test  was  completed  using  the 
bending  channel  only.  The  malfunction  was  in  a  wire  sol¬ 
dered  to  the  terminal,  which  was  repaired  and  caused  no 

further  problems.  A  bending  strain  ga  e  on  the  (-30, /0] 

4  s 
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test  specimen  failed  at  the  end  of  Test  3,  also  due  to  ex¬ 
cessive  deflections.  This  gauge  was  not  repairable,  so 
Tests  14  and  15  were  run  with  only  one  bending  strain  gauge. 
Temperature  fluctuations  during  these  two  tests  caused  a 
zero  drift  on  the  bending  channel.  This  zero  drift  made  the 
data  worthless  for  determining  static  deflections.  However, 
the  data  were  satisfactory  for  determining  oscillation 
frequencies  and  amplitudes. 

The  flutter  threshold  velocity  was  plotted  versus  ini¬ 
tial  angle  of  attack  in  Fig.  8.  The  test  specimens  exhibited 
markedly  different  behavior.  The  [02/90]  test  specimen,  with 
no  bending-torsion  coupling,  exhibited  bending-torsion  flutter 
with  a  frequency  of  approximately  30  Hz  for  angles  of  attack 
below  two  degrees  and  torsion  flutter  with  a  frequency  of 
approximately  42  Hz  for  higher  angles  of  attack.  The  test 
specimens  with  large  positive  bending-torsion  coupling, 

[-*-30 2/°  1  3  anc*  (  +  452/0]s,  exhibited  primarily  bending-torsion 
flutter  at  28  Hz  and  24  Hz,  respectively.  The  flutter  thres¬ 
hold  velocity  did  not  drop  significantly  with  increasing 
angle  of  attack,  probably  because  the  coupling  caused  a  de¬ 
crease  in  the  tip  angle  of  attack,  preventing  it  from  stall¬ 
ing.  At  an  angle  of  attack  of  18  degrees,  the  flutter  thres¬ 
hold  for  the  [+30-/0]  test  specimen  dropped  significantly 
and  the  flutter  changed  to  torsion  flutter  at  55  Hz.  The 
(±45  2/0]  test  specimen,  which  was  stiff  in  torsion  and  had 
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Figure  8.  Flutter  Boundary  Velocities 
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only  a  small  amount  of  positive  bending- torsion  coupling, 
would  not  flutter  at  zero  angle  of  attack  within  the  32  m/sec 
maximum  speed  for  the  wind  tunnel.  However,  at  angles  of 
attack  above  eight  degrees,  it  behaved  similarly  to  the 
[O2/90]s  test  specimen.  The  flutter  was  primarily  torsion 
at  60  Hz.  The  test  specimens  with  large  negative  bending- 
torsion  coupling,  [-3O2/O ]  and  t “ 4 5 2/0 1 s ,  exhibited  divergence 
at  zero  angle  of  attack.  These  two  test  specimens  exhibited 
primarily  bending  flutter  with  frequencies  of  4.5  Hz  and 
4.0  Hz,  respectively,  at  angles  of  attack  greater  than  one 
degree.  This  flutter  was  more  of  a  Von  Karman  vortex  shed¬ 
ding  buffet  then  a  true  stall  flutter. 

Several  data  points  in  Fig.  8  have  been  marked  with  a 
(')  to  indicate  a  high  frequency,  low  amplitude  oscillation. 
The  [O2/90]s  test  specimen  exhibited  bending  oscillations 
at  69  Hz  and  torsion  oscillations  at  145  Hz.  The  [±45_/0] 

4  S 

test  specimen  exhibited  torsion  oscillations  at  108  Hz. 

These  oscillations  indicated  that  the  second  bending  and 
second  torsion  modes  had  been  excited,  which  happened  because 
the  plates  were  quite  flexible.  They  do  not  represent  the 
phenomenon  with  which  this  study  is  concerned. 

The  effect  of  bending-torsion  coupling  on  the  zero 
angle  of  attack  flutter  velocity  could  not  be  determined 
directly  from  Fig.  8.  This  was  because  the  test  specimens 
had  varying  amounts  of  torsional  stiffness,  which  strongly 


affect  the  flutter  velocities.  To  cancel  out  the  effect  of 
torsional  stiffness,  the  dimensionless  reduced  flutter 
velocity  (U_/ba)  )  was  plotted  versus  angle  of  attack  in 
Fig.  9.  In  this  parameter,  b  is  the  semichord  and  is 
the  first  torsion  natural  frequency,  which  was  obtained 
experimentally. 

The  trends  for  each  test  specimen  are  similar  to  Fig. 

8,  as  one  would  expect.  However,  in  Fig.  9,  the  [0-/90] 

fc  s 

test  specimen  had  the  highest  zero  angle  of  attack  reduced 

flutter  velocity,  followed  by  the  [+452/0]s  test  specimen, 

and  finally  the  [+30-/0]  test  specimen.  Since  the 

^  s 

[+302/0]g  test  specimen  had  the  largest  positive  bending- 
torsion  coupling,  the  [+452/0]g  test  specimen  had  somewhat 
less,  and  the  [O2/90]g  test  specimen  had  none,  one  may 
conclude  from  this  very  limited  sample  size  that  zero  angle 
of  attack  flutter  reduced  velocity  was  inversely  propor¬ 
tional  to  positive  bending-torsion  coupling. 

The  [  —  45 2/0 ] s  test  specimen  had  a  higher  reduced  di¬ 
vergence  velocity  (UD/bua)  than  the  [  — 302/0 ]  test  specimen. 
So,  using  the  same  logic  as  before,  one  may  conclude  that 
divergence  velocity  was  inversely  proportional  to  increas¬ 
ingly  negative  bending-torsion  coupling.  Visual  observa¬ 
tions  indicated  that  divergence  was  impossible  for  test 
specimens  with  substantial  positive  bending-torsion 
coupling. 


Reduced  Velocity  ^  U/bw  (dimensionless) 


Figure  9.  Reduced  Flutter  Boundary  Velocities 
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Static  tip  deflections,  both  a. .  and  w  ,  for  the 

fclP  etip 

[+45-/0].  [*45/0]  ,  and  [-45-/0]  were  plotted  versus  velo- 
4  s  s  ^  s 

city  in  Figs.  10,  11  and  12,  respectively.  The  tip  deflec¬ 
tion  data  were  obtained  using  the  data  reduction  procedure 
in  Appendix  E.  This  data  reduction  procedure,  which  relates 
strain  gauge  readings  to  tip  deflections  through  static  cali¬ 
brations,  was  only  an  approximation.  Therefore,  Figs.  10  - 
12  must  be  examined  for  qualitative  information  only.  The 
lines  faired  on  the  plots  represent  a  "best  fit"  to  the 
available  data,  which  have  been  arbitrarily  extrapolated  to 
indicate  some  trends.  The  ot.  .  was  the  sum  of  the  initial 

tip 

angle  of  attack  (a-)  and  the  elastic  twist  (a  ) . 

0  e 

Figure  10  indicates  a  linear  relationship  between 
lateral  tip  deflection  and  velocity.  This  seems  odd,  be¬ 
cause  the  force  on  the  wing  is  a  function  of  the  dynamic 
pressure  (q) ,  which  is  a  function  of  the  square  of  the 
velocity.  One  would  expect  this  relationship  to  be  para¬ 
bolic.  Two  possible  explanations  are  offered.  First,  the 
relationship  was  parabolic  and  the  data  were  somewhat  in 
error.  Second,  the  bending- torsion  coupling  reduced  the 
tip  angle  of  attack  to  a  low  value  as  the  test  specimen  was 
bent  laterally.  Unfortunately,  the  plot  of  the  tip  angles 
of  attack  doesn't  support  this  second  explanation.  However, 
values  for  obtained  from  the  strip  chart  were  extremely 


Figure  10.  Experimental  Static  Aeroelastic  Deflection,  [+45 VO]  Te 
Specimen  s 
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Figure  11.  Experirrental  Static  Aeroelastic  Deflections,  [±45/0]s  Itest 
Specimen 
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sensitive  to  slight  errors  in  reading  the  strip  chart.  This 
was  compounded  by  the  fact  that  buffeting  of  the  test  speci¬ 
mens  made  it  difficult  to  get  an  accurate  strip  chart  reading. 

Figure  11  indicates  a  more  parabolic  relationship  be¬ 
tween  the  lateral  tip  deflections  and  velocity.  Also  the  tip 
angle  of  attack  increases  with  increasing  velocity  as  one 
would  expect  for  a  wing  with  low  bending-torsion  coupling. 

More  data  points  were  available  for  Fig.  11  and  they  have 
good  continuity. 

Figure  12  shows  the  expected  divergent  response  ten¬ 
dencies  of  a  wing  with  negative  bending- torsion  coupling. 

The  flutter  amplitude  was  plotted  versus  velocity  for 
the  [+45-/01  ,  [*45/0]  ,  and  [-45-/0]  test  specimens  in 

S  S  c,  S 

Figs.  13,  14,  and  15,  respectively.  The  flutter  amplitude 
was  unreduced,  unamplified,  peak  to  peak  oscillation  data 
directly  from  the  strip  chart.  It  has  no  significance 
other  than  to  indicate  where  the  flutter  starts,  for  a  given 
initial  angle  of  attack,  and  to  indicate  how  rapidly  the 
magnitude  increases.  For  both  the  [  +  452/0Js  and  £  —  452/0 ] s 
test  specimens  flutter  beings  quite  close  together  for 
several  angles  of  attack.  While,  for  the  [±45/0]g  test 
specimen  the  occurrence  is  more  spread  out. 


Strip  Chart  Peak-to-Peak  (Bending)  Anplitude  ^  A  (nm) 


Strip  Chart  Peak-to-Peak  (Torsicn)  Anplitude  %  A  (inn) 
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CHAPTER  IV 

COMPARISON  OF  THEORETICAL  AND  EXPERIMENTAL  RESULTS 

1 8  19 

Work  done  by  Boyce  ,  Crawley  ,  and,  in  particular, 
Turner10  showed  that  stiffness  properties  for  graphite/epoxy 
laminates  experimentally  obtained  from  in-plane  extensional 
loading  were  generally  higher  than  material  properties  ob¬ 
tained  from  out-of-plane  flexural  loading.  The  standard 
TELAC  values  for  E^,  E^,  GLT  and  v^,  which  appear  in  Ap¬ 
pendix  A,  were  obtained  for  Hercules  AS/3501-6  G/E  from  in¬ 
plane  loads.  Turner,  through  a  series  of  out-of-plane 
flexural  static  and  dynamic  tests,  established  an  alterna¬ 
tive  set  of  material  property  values  for  AS/3501-6  G/E. 

These  material  property  values  also  appear  in  Appendix  A. 
Turner’s  values,  while  not  statistically  significant  like 
the  TELAC  values,  do  constitute  a  basis  for  comparison. 
Therefore,  in  this  chapter,  experimental  results  will  be 
compared  to  theoretical  results  using  both  TELAC  and  Turner 
material  properties  for  graphite/epoxy.  Additionally,  dis¬ 
cussion,  in  several  instances,  will  be  limited  to  quali¬ 
tative  "trends"  due  primarily  to  the  approximate  nature  of 
the  data  being  compared. 

4.1  Static  Deflection  Tests 


The  theoretical  and  experimental  flexibility  influence 
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coefficients  for  each  test  specimen  appear  in  Table  3.  For 
clarity,  they  have  been  arranged  in  2  x  2  matrices.  As 
stated  before,  c^  is  the  bending  flexibility,  c^2  is  the 
torsion  flexibility,  and  c^2,  are  the  bending  torsion 

coupling  flexibilities. 

Values  for  c^1  calculated  from  Turner  material  proper¬ 
ties  were  within  10%  of  experimental  values.  The  c^'s 
from  TELAC  material  properties  were  lower  than  Turner 
and  experimental  c-^'s,  indicating  that  TELAC  material  pro¬ 
perties  resulted  in  a  plate  that  was  too  stiff  in  pure 
bending.  The  coupling  flexibility  influence  coefficients 
obtained  from  TELAC  data  were,  however,  in  better  agreement 
with  experimental  data  than  those  from  Turner  data. 

The  torsion  influence  coefficient  for  the  [*45/0]  test 

s 

specimen  obtained  from  Turner  data  had  reasonable  agreement 
with  the  experimental  value.  The  torsion  influence  coef¬ 
ficients  obtained  from  both  Turner  and  TELAC  data  had  poor 
agreement  with  experimental  values  for  test  specimens  with 
large  bending- tors ion  coupling.  With  one  notable  exception, 
all  theoretical  values  for  c22  were  lower  than  experimental 
values.  The  exception  was  the  [O2/90 ]  test  specimen,  where 
theory  predicted  a  substantially  higher  value  for  c22  than 
was  obtained  from  the  experiments.  For  all  test  specimens, 
except  [O2/90]s,  torsion  influence  coefficients  obtained 
using  Turner  data  had  better  agreement  with  experimental 
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values  than  did  torsion  influence  coefficients  obtained 
us ing  TELAC  data . 

There  were  several  possible  reasons  for  differences 
between  theoretical  and  experimental  torsion  influence  co¬ 
efficients.  First,  the  two  term  deflection  equation  may 
not  adequately  represent  the  actual  deflection  of  the 
structure.  Additional  terms  would  allow  more  degrees  of 
freedom.  Second,  the  assumption  that  the  plate  was  chord- 
wise  rigid  when  undergoing  torsional  deflections  may  not 
be  true,  especially  for  test  specimens  with  large  bending- 
torsion  coupling.  Third,  the  theory  used  to  obtain  the 
flexural  moduli  (D.^)  for  the  plates,  as  presented  in 
Chapter  II,  assumed  small  deflections  and,  consequently, 
neglected  transverse  shear.  If  transverse  shear  was  pre¬ 
sent,  the  plate  would  become  more  flexible  in  torsion. 
Further,  if  transverse  shear  increased  with  increasing 
bending- torsion  coupling,  then  for  test  specimens  with 
large  bending-torsion  coupling,  one  would  expect  the  experi¬ 
mental  torsion  flexibility  influence  coefficients  to  be 
much  larger  than  the  theoretical  influence  coefficients. 

This  was  the  case  in  this  study. 

The  final  significant  trend  was  that  the  coupling 
factor  (k)  for  the  theoretical  flexibility  influence  co¬ 
efficients  was  substantially  higher  than  the  experimental 
coupling  factor  for  all  test  specimens  except  [Oj/SOlg. 
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Both  the  theoretical  and  experimental  coupling  factors  were 
essentially  zero  for  the  [C>2/90]s  test  specimen. 


4.2  Free  Vibration  Tests 


A  summary  of  the  theoretical  and  experimental  natural 

frequencies  for  the  first  three  natural  vibration  inodes  of 

each  test  specimen  appears  in  Table  4.  Also  included  in 

Table  4  are  the  theoretical  uncoupled  natural  frequencies 
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calculated  using  the  following  equations. 


(first  bending) 


(second  bending) 


(first  torsion) 


“IB 


“2B 


“IT 


1.57 


48D66 

m£2c2 


(4-la) 

(4-lb) 

( 4-lc) 


where  is  the  bending  flexural  modulus,  Dgg  is  the  tor¬ 
sion  flexural  modulus,  m  is  the  mass  per  area,  l  is  the 
length,  c  is  the  chord,  and  oj  is  the  frequency  in  rad/sec. 

The  theoretical  first  bending  natural  frequencies, 
calculated  using  Turner  material  properties ,  were  all  within 
5%  of  the  experimental  values,  which  was  good  agreement. 

The  first  bending  frequencies  calculated  from  TELAC  material 
properties  had  poorer  agreement  with  experimental  values. 
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In  all  cases,  the  TELAC  first  bending  frequencies  were 
higher  than  experimental  frequencies.  The  theoretical  un¬ 
coupled  first  bending  frequencies,  as  one  would  expect, 
showed  good  correlation  with  experimental  frequencies  for 
test  specimens  with  little  bending- torsion  coupling,  and 
poor  correlation  for  test  specimens  with  large  bending- 
torsion  coupling. 

The  theoretical  values  for  the  second  bending  fre¬ 
quencies  from  Turner  material  properties  showed  reasonable 
agreement  with  the  experimental  values  for  test  specimens 
with  low  coupling,  and  poor  agreement  for  test  specimens 
with  large  coupling.  Second  bending  frequencies  calculated 
from  TELAC  material  properties  were  not  as  good  as  the 
frequencies  from  Turner  data.  The  theoretical  uncoupled 
values  for  the  second  bending  frequency  had  generally 
poorer  agreement  with  the  experimental  values  than  the 
coupled  solutions  did. 

The  theoretical  first  torsion  natural  frequencies  had 
generally  poor  agreement  with  experimental  values.  Errors 
for  first  torsion  natural  frequencies  calculated  using 
Turner  material  properties  were  lowest  for  the  [ *4 5/0 ] g  test 
specimen  and  highest  for  the  [+452/0]s  test  specimen.  More 
importantly,  one  observes  that  the  theoretical  torsion  fre¬ 
quency  values  for  the  [^45/0]  and  [+45-/0]  test  specimens 
were  essentially  the  same,  while  the  experimental  torsion 


frequency  for  the  t+452/0]g  test  specimens  was  substantially 
lower  than  the  experimental  torsion  frequency  for  the  [*45/0] 
test  specimen.  The  higher  torsion  frequency  for  the  [*45/0 ] 
test  specimen  indicated  that  it  was  much  stiffer  in  torsion 
than  the  [+45-/0]  test  specimen.  As  before,  this  was  pos- 

Z  S 

sibly  due  to  neglecting  transverse  shear  in  the  problem 
derivation.  Also,  the  fact  that  only  the  first  three  vib¬ 
ration  mode  shapes  were  used  to  formulate  the  assumed  deflec¬ 
tion  equation  could  affect  the  solutions  somewhat.  Finally, 
as  before,  the  assumption  that  the  plate  was  chordwise  rigid 
during  torsion  vibrations  could  also  affect  the  solution. 

The  theoretical  natural  frequencies  should  possibly  be 
checked  by  a  more  extensive  Rayleigh-Ri tz  formulation  and 
a  finite  element  analysis. 

4 . 3  Divergence  Velocities 

Theoretical  and  experimental  divergence  velocities  for 
each  test  specimen  appear  in  Table  5.  Velocities  were  ob¬ 
tained  for  plates  using  both  TELAC  and  Turner  material  prop¬ 
erties,  and  for  both  a  two  dimensional  lift  curve  slope 
(2c^/2a)  and  one  empirically  corrected  for  finite  span  using 
the  correction  expression  in  Chapter  II.  Of  the  six  speci¬ 
mens  tested,  only  three  would  possibly  exhibit  any  divergence 
at  zero  initial  angle  of  attack.  These  three  test  specimens 
were  [O2/90]s,  [  —  45 2/0 ] g ,  and  [-302/0]g. 


Of  these  three. 
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the  [0-/90]  test  specimen  exhibited  flutter  prior  to  diver- 

^  S 

gence,  so  no  experimental  divergence  characteristics  were 
obtained.  However,  a  theoretical  divergence  velocity  was 
still  calculated  for  that  specimen. 

Theoretical  divergence  velocities  were  lower  than  experi¬ 
mental  velocities  in  all  cases.  The  finite  span  divergence 
for  the  [-452/0]s  and  [-302/0]g  test  specimens,  calculated 
using  Turner  material  properties,  were  lower  than  the  diver¬ 
gence  velocities  calculated  using  TELAC  material  properties. 

An  experimental  divergence  velocity  was  not  obtained  for 
the  [O2/90]g  test  specimen.  However,  the  wind  tunnel  was 
taken  up  to  27  m/sec  (90  ft/sec)  with  this  test  specimen 
and  divergence  did  not  occur.  This  contrasts  to  a  theoreti¬ 
cal  finite  span  divergence  velocity,  calculated  with  Turner 
material  properties,  of  22  m/sec  (71  ft/sec) .  This  indi¬ 
cated  that  the  theory  would  not  adequately  model  the  [O2/90]g 
test  specimen. 

The  fact  that  the  theory  predicted  a  lower  divergence 
velocity  than  found  from  experiments  would  suggest  that  the 
actual  plates  were  stiffer  than  the  theoretical  plates.  This 
observation  is  in  sharp  contrast  to  the  observations  in  Sec¬ 
tions  4.1  and  4.2,  where  the  theoretical  plates  appeared 
stiffer.  The  observation  would  tend  to  rebut  the  explana¬ 
tions  given  in  those  sections.  Since  divergence  is  a  func¬ 
tion  of  aerodynamic  force  on  the  plate,  as  well  as  plate 


9? 


stiffness,  the  aerodynamic  strip  theory  needed  to  calculate 
the  loads  on  the  plate  may  have  been  a  source  of  error. 
Further,  the  finite  span  correction  may  have  resulted  in  the 
theoretical  air  loads  being  greater  than  the  actual  air  loads. 
Reference  13  suggests  that  an  alternate  empirical  correction 
for  finite  span  effects  may  provide  better  correlation  be¬ 
tween  theoretical  and  experimental  results  for  divergence 
problemr.  That  empirical  correction  is 

— ^  a  a-  E  — — —  ]  (4-2) 

3a  AR  +  4 

Using  this  finite  span  correction,  theoretical  divergence 
velocities  were  raised  slightly.  A  final  source  of  error 
may  have  been  the  criteria  used  during  the  wind  tunnel 
tests  to  determine  the  divergence  velocity.  This,  however, 
is  an  improbable  source  of  error. 

4 . 4  Static  Aeroelastic  Tip  Deflections 

The  theoretical  tip  lateral  deflections  (w-,  )  and 

t-Tip 

tip  angles  of  attack  (aTIp)  were  calculated  using  Turner 
material  properties,  the  empirical  finite  span  correction, 
and  the  solution  technique  in  Chapter  II.  The  tip  deflec¬ 
tions  were  plotted  versus  velocity  for  the  [+45-/0] 

S 

[-45/0]g,  and  [-452/0]g  test  specimens  in  Figs.  16,  17,  and 
18,  respectively.  To  allow  a  comparison  with  experimental 
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results,  the  approximate  experimental  tip  deflections,  dis¬ 
cussed  in  Chapter  III,  were  included  in  these  figures.  The 
lines  represent  the  theory  and  the  plotting  symbols  are  the 
experimental  values.  Since  data  from  both  theory  and  experi¬ 
ment  were  approximate,  a  valid  comparison  must  be  limited 
to  trends  only. 

In  Fig.  16,  the  lateral  deflection  trends  for  the 

[+45-/0]  test  specimen  show  reasonable  agreement  between 
z  s 

theory  and  experiment,  especially  at  higher  velocities. 

The  trends  for  the  theoretical  and  experimental  tip  angle 
of  attack,  however,  did  not  agree.  The  experimental  values 
indicated  a  very  slight  increase  in  tip  angle  of  attack  as 
velocity  was  increased.  Theory,  on  the  other  hand,  showed 
a  substantial  decrease  in  tip  angle  of  attack  as  velocity 
was  increased.  The  implications  of  the  theoretical  results 
were  that  the  torsional  moment  resulting  from  bending-torsion 
coupling  more  than  overcame  the  aerodynamic  moment.  A  pos¬ 
sible  reason  was  that  the  theoretical  coupling  factor  (k) , 
discussed  in  Section  4,1,  was  much  larger  than  the  experi¬ 
mental  k .  Of  course,  there  was  also  the  problem  with  read¬ 
ing  the  strip  chart,  as  discussed  in  Chapter  III. 

The  agreement  between  theoretical  and  experimental 
lateral  tip  deflections  for  the  [*45/0]  test  specimen  in 

S 

Fig,  17  was  not  quite  as  good  as  for  the  [  +-4 5 2/°  1  s  test 
specimen,  but  was  still  reasonable.  Again,  the  trends  for 


1 


2 


the  tip  angles  of  attack  did  not  correlate  between  theory 

and  experiment.  The  theoretical  tip  angles  of  attack  did 

not  decrease  as  much  for  the  [*45/0]  test  specimen  as  they 

did  for  the  [+45-/0]  test  specimen,  but  the  experimental 

fc  s 

tip  angle  of  attack  values  increased  more  for  the  [*45/0] 

s 

test  specimen.  In  view  of  the  fact  that  the  [*45/0 ]  test 
specimen  possessed  slight  positive  bending- torsion  coupling, 
one  would  anticipate  this  type  of  behavior.  The  reasons 
for  the  lack  of  correlation  between  experiment  and  theory 
then  would  be  the  same  as  for  the  [+452/0 ]g  test  specimen. 

Divergent  tendencies  of  both  the  tip  lateral  deflec¬ 
tions  and  angles  of  attack  for  the  [-452/0]g  tesh  specimen 
are  apparent  in  Fig.  18.  The  trends  were  similar  for  both 
the  theoretical  and  experimental  results,  although  the  ex¬ 
perimental  data  did  not  correlate  well  with  theory  for  indi¬ 
vidual  angles  of  attack.  Finally,  the  experimental  tip  angle 
of  attack  data  appeared  to  diverge  faster  than  the  theore¬ 
tical  data. 

4 . 5  Flutter  Velocities 

The  theoretical  and  experimental  flutter  velocities  at 
zero  angle  of  attack  for  each  test  specimen  appear  in  Table 
6.  Additionally,  the  nondimensional  reduced  flutter  velo¬ 
cities  for  each  test  specimen  are  presented.  Since  the 
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TABLE  6 

a.  FLUTTER  VELOCITIES 


Test 

Specimen 

Theoretical 

TELAC 

(U-g  Analysis) 
Turne  r 

Experimental 

" 

[o2/9°]s 

18,0 

(59,0 

m/sec 
f t/secl 

17,4  m/sec 
(57.0  ft/sec) 

25.9  m/sec 
(85.0  ft/sec) 

[±45/0] s 

42.1 

(138) 

36.6 

(120) 

32.0 

(>105) 

[+452/°]s 

40.2 

(132) 

35.1 

(115) 

27.7 

(91.0) 

t-452/0]s 

50.3 

(165) 

43.6 

(143) 

diverged 

first 

[+302/0]s 

34.7 

(114) 

30.5 

(100) 

27.1 

(89.0) 

[-303/0]s 

50,0 

(164) 

42.3 

(142) 

diverged 

first 

b.  REDUCED  FLUTTER  VELOCITIES 

Test 

Specimen 

TELAC 

Theoretical 

Turner 

Experimental 

[o2/9°]s 

2.21 

2.21 

2.56 

[*45/0 ]g 

2.21 

2.21 

>1.6 

[+452/°]s 

2.11 

2.11 

2.24 

[-45^/0] 
l  s 

2.64 

2.63 

diverged 
first  (1.01) 

[+302/01s 

2.04 

2.04 

1.94 

C-302/0]a 


2.93 


2.90 


diverged 
first  (.839) 
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[-452/O]  and  I-302/0Js  test  specimens  diverged  prior  to 
fluttering,  no  experimental  flutter  velocities  were  obtained 
for  them.  However,  theoretical  flutter  velocities  were 
calculated  for  these  two  test  specimens.  Also,  the  [±45/0]s 
test  specimen  did  not  flutter  within  the  speed  range  of  the 
wind  tunnel. 

The  correlation  between  theoretical  and  experimental 
flutter  velocities  was  generally  poor.  The  flutter  velocity 
for  the  [O2/90]g  test  specimen,  calculated  using  Turner 
material  properties,  was  much  lower  than  the  experimental 
flutter  velocity.  On  the  other  hand,  the  theoretical  flutter 
velocities  for  the  [+452/0]g  and  [+302/0]g  test  specimens, 
calculated  using  Turner  material  properties,  were  much 
higher  than  the  experimental  flutter  velocities.  The  flutter 
velocities  calculated  using  TELAC  material  properties  had 
generally  worse  correlation. 

The  errors  in  the  flutter  calculations  appeared  to  be 
directly  related  to  the  inability  of  the  theory  to  accurately 
determine  the  torsion  vibration  frequencies.  Therefore,  one 
would  expect  a  somewhat  better  correlation  between  theoreti¬ 
cal  and  experimental  reduced  flutter  velocities  (Up/bw^) , 
since  the  contribution  of  the  torsion  vibration  frequency 
had  been  eliminated.  In  fact,  that  was  exactly  what  hap¬ 
pened.  The  [02/90]s  test  specimen's  error  was  reduced  sig¬ 
nificantly,  and  the  errors  associated  with  the  [+45-/0] 

2  s 
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and  [+30-/0]  test  specimens  were  less  than  10%.  The 
z  s 

reader  should  note  that  the  ui  used  for  the  theoretical 
calculations  was  obtained  from  the  U-g  analsyis  for  a  very 
large  reduced  frequency  (k)  value. 

Since  the  trends  were  reasonably  well  predicted  by  the 
theory,  some  conclusions  on  how  bending- torsion  coupling 
affects  flutter  velocity  were  drawn,  based  on  theoretical 
results.  The  first  conclusion  was  that  increasingly  posi¬ 
tive  bending-torsion  coupling  reduced  the  flutter  velocity. 
The  second  conclusion  was  that  increasingly  negative  bending- 
torsion  coupling  raised  the  flutter  velocity.  This  was  an 
unrealistic  situation,  however,  because  the  test  specimens 
with  negative  coupling  diverged  long  before  they  reached 
the  flutter  velocity. 
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CHAPTER  V 

CONCLUSIONS  AND  RECOMMENDATIONS 

5 .1  Conclusions 

The  present  study  has  sought  to  investigate  the  flutter 
and  divergence  behavior  of  graphite/epoxy  plates  with  various 
amounts  of  bending-torsion  stiffness  coupling.  As  a  conse¬ 
quence,  the  static  deflection  and  vibration  behavior  has  also 
been  examined. 

There  was  reasonable  agreement  between  theory  and  exp¬ 
eriments  for  the  bending  stiffness  and  first  bending  mode 
natural  frequencies.  Also,  for  test  specimens  with  low 
bending- torsion  coupling,  there  was  reasonable  agreement 
for  the  torsional  stiffness  and  first  torsion  mode  natural 
frequencies.  There  was,  however,  poor  agreement  between 
theory  and  experiments  for  torsion  stiffness  and  frequencies 
of  test  specimens  with  large  bending-torsion  coupling.  The 
poor  agreement  was  probably  due  to  an  inadequate  number  of 
modes  in  the  assumed  deflection  equation.  Other  possible 
sources  of  error  were  the  assumption  that  the  plate  was 
chordwise  rigid  and  the  neglecting  of  the  transverse  shear 
effect. 

In  this  study,  only  the  plates  with  negative  bending- 
torsion  coupling  (-D^g)  exhibited  divergence  within  the 
speed  range  of  the  wind  tunnel.  For  these  test  specimens, 
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the  theory  had  reasonable  agreement  with  the  experiments. 

The  errors  in  the  theory  were  possibly  due  to  the  approxi¬ 
mate  aerodynamic  strip  theory.  For  test  specimens  with 
positive  bending- torsion  coupling,  the  theory  predicted  that 
divergence  would  not  occur.  Experimentally,  they  did  not 
exhibit  a  divergence  tendency  within  the  speed  range  of  the 
wind  tunnel.  Based  on  the  theory  and  the  limited  set  of 
experiments,  the  following  observations  can  be  made.  The 
zero  angle  of  attack  divergence  velocity,  for  the  plates  in 
this  study,  became  infinite  for  positive  bending-torsion 
coupling  above  a  certain  critical  value.  This  divergence 
velocity  decreased  as  the  bending-torsion  coupling  decreased 
from  the  critical  value  to  zero,  and  continued  to  decrease 
as  the  bending- torsion  coupling  became  increasingly  negative. 

The  theory  and  experiments  had  poor  agreement  for  zero 
angle  of  attack  flutter  velocities.  However,  the  agreement 
between  theory  and  experiments  for  reduced  flutter  velo¬ 
cities  (U/bw^)  was  good.  This  indicated  that  the  inaccu¬ 
racies  in  the  theory  were  in  the  prediction  of  the  torsion 
vibration  frequency.  The  test  specimens  with  negative 
bending-torsion  coupling  diverged  prior  to  reaching  their 
flutter  velocities,  but  a  theoretical  flutter  velocity  was 
still  calculated  for  them.  As  before,  based  on  the  theory 
and  the  limited  set  of  experiments,  the  following  observa¬ 
tions  can  be  made . 
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The  zero  angle  of  attack  flutter  velocity  for  the  plate  in 
this  study  was  low  for  positive  bending-torsion  coupling. 

The  flutter  velocity  increased  as  the  positive  bending- 
torsion  coupling  decreased  to  zero,  and  continued  to  in¬ 
crease  as  the  bending-torsion  coupling  became  increasingly 
negative.  This  latter  situation,  for  negative  bending- 
torsion  coupling,  becomes  academic  if  divergence  occurs 
prior  to  flutter. 

The  test  specimens  exhibited  markedly  different  flutter 
behavior  as  angle  of  attack  was  increased.  The  test  speci¬ 
mens  with  low  bending-torsion  coupling  exhibited  classical 
bending- torsion  flutter  at  zero  angle  of  attack,  but  it 
changed  to  pure  torsion  flutter  as  the  angle  of  attack  was 
increased.  Also,  the  flutter  boundary  velocity  decreased 
rapidly  with  increasing  angle  of  attack,  stabilizing  around 

g 

a  =  12°,  This  trend  is  similar  to  the  results  Rainey  ob¬ 
tained  for  metal  plates  and  is  classical  stall  flutter.  The 
test  specimens  with  large  positive  bending-torsion  coupling 
exhibited  bending-torsion  flutter  throughout  most  of  the 
angle  of  attack  range.  Additionally,  the  flutter  boundary 
velocities  did  not  decrease  with  increasing  angle  of  attack 
below  a  =  12°.  This  behavior  was  probably  due  to  the  fact 
that  the  tip  of  the  plate  tended  to  have  a  lower  angle  of 
attack  than  the  root,  so  it  did  not  stall.  The  test  specimens 
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with  large  negative  bending-torsion  coupling  diverged  at 
zero  angle  of  attack  and  exhibited  a  low  frequency  bending 
flutter  for  angles  of  attack  greater  than  zero.  The  flutter 
boundary  velocities  here  were  significantly  below  the  velo¬ 
cities  for  the  other  test  specimens  for  angles  of  attack 
less  than  12°.  These  bending  flutter  velocities  decreased 
moderately  with  increasing  angle  of  attack. 

5 . 2  Recommendations 

The  theory  should  be  further  developed  to  give  better 
results  for  the  torsion  vibration  frequencies  and  flutter 
velocities.  Improvement  could  probably  be  most  easily  faci¬ 
litated  by  going  to  a  finite  element  analysis  which  allows 
more  degrees  of  freedom  and  accounts  for  transverse  shear. 

If  a  Rayleigh-Ritz  formulation  is  used,  more  terms  should 
be  included  in  the  deflection  equation.  For  divergence 
velocity  calculations,  a  better  aerodynamic  theory,  such 
as  lifting  line  theory,  should  be  used.  Also,  the  stall 
characteristics  of  the  plate  should  be  explored  further. 

The  effect  of  bending-torsion  coupling  as  a  passive 
gust  alleviation  technique  might  hold  some  promise.  Further, 
bending-torsion  coupling  might  be  useful  for  increasing  the 
aileron  reversal  speed  for  a  swept  wing.  Finally,  the  effect 
of  bending-torsion  coupling  on  an  oblique  wing,  particularly 
with  the  roll  degree  of  freedom,  might  prove  interesting. 
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APPENDIX  A 

MATERIAL  PROPERTIES 


TABLE  A-l 

ORTHOTROPIC  ENGINEERING  CONSTANTS 
Hercules  AS/3501-6  Graphite/Epoxy 


Property 


TELAC  Values  Turner  Values 

Cin-plane  loads)  (out-of-plane  loads) 


°LT 


Ply  thickness 
Density 


130  x  109  Pa 
10.5  x  109 
0.28 

6.00  x  109 

1.34  x  10-6  m 
1.52  kg/m9 


98  x  109  Pa 
7.9  x  109 
0.28 

5.6  x  109 


Ill 


TABLE  A- 2 

LAMINATE  FLEXURAL  MODULI 


Laminate  Modulus  TELAC  Turner 

Material  Properties  Material  Properties 


D11 

5.474  Nm 

“  ‘  .  1 

4.038  lb • ft 

4.126  Nm 

3.043  lb • ft 

[02/9°]s 

°16 

0 

0 

0 

0 

°66 

0.2600 

0.1918 

0.2425 

0.1789 

°11 

1,996 

1.472 

1.550 

1.143 

[*45/0 ] 

s 

°16 

0.5789 

0.427 

0.4364 

0.3219 

°66 

1.422 

1.049 

1.074 

0.7920 

°11 

1.996 

1.472 

1.550 

1.143 

t+4V0]s 

°16 

1.254 

0.925 

0.9456 

0.6975 

°66 

1.422 

1.049 

1.074 

0.7920 

D11 

3.541 

2.612 

2.703 

1.994 

[+302/0]s 

°16 

1,589 

1.172 

1.179 

0.8695 

D66 

1.132 

0.8350 

0.8662 

0.6389 
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B.l  Mode 
*x(x)  = 

where  £^  = 
al 

<p2(x)  = 

where  e2 
a2 

<t>1(Z)  = 

<P2  ( z)  = 

»a(*>  - 


APPENDIX  B 
THEORETICAL  ANALYSIS 


Shapes 


14,15 


cosh£l(j)  -  cos£;l(j)  -  ai  [sinh£l  (2)  -  sine1(j)] 

(B-l) 


1.8751 

0,7341 


coshe2(^)  -  cose2(j)  -  a2Isinhe2{f)  “  sine2(f)l 

(B-2) 

4.6941 

1.0185 

<J>  Cx)  =  sin(^y)  (B-3) 

21 


2.00 

4>1  (0 . 75fi.) 

=  1.32 

2.00 

<j>2  (0 . 75 £,) 

=  0.269 

.00 

<P  (0.75  Z) 
ra 

=  0.924 

The  following  integrals  were  evaluated  numerically  using 
ten  point  trapezoidal  integration. 


1  rl  *  a 

7  J  4>-,dx 

*  0  1 


1  r*  *  j 

1  /  *adx 

0 


=  0.783 


=  0.637 


J3  =  i  /  *1*«  dx  =  °-678 


I4  H  I  0/A  *1  dx  =  1*00 


I  ^  *a  dx  =  0*500 


£2  of  *£*a  dx  *  3’76 


£3  g7  (<t>”)2  dx  =  12.4 


1  0f  (<^)2  dx  =  1.23 


12  Qf  *2*o  dx  =  _6’43 


£3  O^  (<^25  2  dx  =  485 


11  = 


'12  = 


</  (<*>")  dx  =  3,04 


t  </  *2  d*  *  i-00 


:ural  Danping  ^  g  (dimensionless)  Structural  Danping  ^  g  (dimensionless) 


w 


«  . 


Figure  B-l.  U-g  Diagram,  [O^SOJg  Test  Specimen 


Figure  B-2.  U-g  Diagram,  [*45/0]  Test  Specimen 


Structural  Danping  ^  g  (dinensionless)  Structural  Danping  ^  g  (dimensionless) 
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Figure  B-4.  U-g  biagram,  [+30^/0 ]g  Test  Specimen 


f — r 


Structural  Danping  ^  g  (dimensionless)  Structural  Danping  ^  o  (dimensionless) 


APPENDIX  C 

GRAPHITE/EPOXY  CURING  TECHNIQUES 


3  Vacuum  bag 

Air  bleeder 
a  (fiberglass  cloth) 

3  Porous  Teflon 

1  Aluminum  plate 

J  (top) 
a  Nonporous  Teflon 

Paper  bleeder 

3  Porous  Teflon 


6-ply  G/E  laminate 
(with  peel  ply) 

Porous  Tteflon 

Paper  bleeder 

Nonporous  Teflon 

Aluminum  plate 
(bottom) 


Figure  C-l.  Curing  Materials 
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TABLE  C-l 

AUTOCLAVE  CURE  CYCLE 


1.  Place  aluminum  curing  plate  with  curing  materials  (Fig. 
C-l)  in  autoclave  and  attach  vacuum  hose. 

2.  Perform  vacuum  tests  Apply  full  vacuum  (30  in)  to  curing 
plate,  turn  off  vacuum  pump,  and  watch  bleed-down  for 
five  minutes. 

3.  Seal  autoclave,  apply  full  vacuum,  and  raise  autoclave 
pressure  to  85  psi. 

4.  Raise  autoclave  temperature  to  240°  F. 

5.  Maintain  85  psi  and  240°  F  for  one  hour. 

6.  Raise  pressure  to  90  psi. 

7.  Raise  temperature  to  35Q°  F. 

8.  Maintain  90  psi  and  350°  F  for  two  hours. 

9.  Decrease  temperature  at  a  rate  not  exceeding  5°  per 

minute  until  temperature  reaches  130°  F. 

10.  Release  pressure  and  vacuum. 

11.  Open  autoclave  and  remove  curing  plate. 
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APPENDIX  D 

TEST  SPECIMEN  DIMENSIONS 

TABLE  D-l 

NOMINAL  TEST  SPECIMEN  DIMENSIONS 


Parameter  Value 

Length  12,0  in  304.8  mm 

Width  3.00  76.2 

Ply  thickness  0.00528  0.134 

(nominal  TELAC  value) 

Total  thickness  0.0317  0.804 

(6  x  ply  thickness) 


Density  0.055  lb/in3  1.52  x  103  kg/m3 

(nominal  TELAC  value) 


EXPERIMENTAL  TEST  SPECIMEN  AVERAGE  MEASUREMENTS 
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APPENDIX  E 

TEST  DATA  AND  DATA  REDUCTION  METHODS 

E. 1  Flexibility  Influence  Coefficient  Data  Reduction 

Test  apparatus  characteristics: 

moment  arm  (1/2  couple)  =  3.5  in 

distance  between  front  and  rear  rulers  =  8.0  in 

Test  1  raw  data: 

applied  forces  (F^) :  ±  0,022  lb,  ±  0.044,  ±  0.066, 

±  Q.088 

obtain  front  and  rear  deflections  for  each  applied  force 
Test  2  raw  data: 

applied  moments  (M^) :  ±  0,013  ft*lb,  ±  0.026,  ±  0.044, 

±  0.056 

obtain  front  and  rear  deflections  for  each  applied  moment 
Data  reduction  for  each  test  point: 

front  +  rear 

lateral  deflection  (w„)  in  units  of  feet  =  [  - 

E  768 

angular  deflection  (a)  in  units  of  radians  = 

front  -  rear 


tan  ^  [ 


256 


] 
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wT 


'11 


'21 


a  WE 

;  C12  *  ;  C22 


a 

Ml 


Plots  of  all  test  points  for  each  test  specimen  appear  in 
Figs.  E-l  to  E-4. 

E. 2  Indicated  and  True  Velocity  Calculations 
Raw  data: 

manometer  reading  (hal(=)  in  inches  of  alcohol 
temperature  (°F)  at  start  and  finish 
barometric  pressure  (h^)  in  inches  of  mercury 

Indicated  velocity  calculation: 

specific  gravity  (Sg)  of  alcohol  =  0.805  at  80°  F 
specific  weight  (u)alc)  Of  alcohol  =  Sg(a)H  Q)  = 

(0,805)  (62.4  lb/ft3)  =  50.2  lb/ft3 


U.  . . 
indicated 


2uialchalc 

PSL 


=  59.3 


halc (inches) 


True  velocity  calculation: 


“true  “  “indicated1  >  "here  ®  * 

PSL  1 


p/pSL  and 
0.002376  sl/ff 


p  =  Pa(psf)/RT  (°K) 

d  a 


i 


J 


M  (ft/lb) 


Figure  E-2.  Flexibility  Influence  Coefficients,  [*45/0]  Test  Specimen 


Figure  E-4.  Flexibility  Influence  Coefficients,  [+302/0]g  Test  Specimen 
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R  »  3089.8  ft/sec2°K 

T  ( °K)  *  0 .5556  [  {  ~  T-St°P  )  +  460  ] 

a  2 

Pa(psf)  =  SgHgu)H^0hHg/12  =  70.72  hHg  (inches) 

where  Sg„  =  13.6 
tig 

Utrue  =  59 • 3  /halc(inches)/a 


E. 3  Experimental  Tip  Deflection  Calibrations 


Calibration  1  raw  data: 

tip  lateral  deflections  (6A) :  ±  0.25  in;  ±  0.5;  ±  0.75; 

±  1.0 

obtain  strain  gauge  readings  from  strip  chart  in  units 
of  mm 

bending  strain  gauge  (channel  2) 
torsion  strain  gauge  (channel  1)  (0c) 


Calibration  2  raw  data: 

tip  angular  deflections  (eA>:  ±  4.78  deg;  ±  9.55 
obtain  strain  gauge  readings  from  strip  chart  in  units 
of  nun 

By  plotting  each  calibration  point,  obtain  the  following 


slopes : 
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all 

a21 

a12 

a22 


5c 

(cal 

1) 

5A 

(cal 

1) 

9c 

(cal 

1) 

5A 

(cal 

1) 

5c 

(cal 

2) 

0A 

(cal 

2) 

9c 

(cal 

2) 

9A 

(cal 

2) 

From  these,  we  can  obtain  the  following  transformation  equa¬ 
tion: 


S (mm) 
c 

^  _ 

"all 

a12 

I  6A(in) 

0  (mm) 
c 

a21 

a22 

j  0A(deg) 

Since  we  want  the  inverse  relationship,  the  A-matrix  is 
inverted  to  obtain  the  final  expression. 


«  [  A  ] 


-1 


From  this  expression,  we  can  get  approximate  values  for  the 
test  specimen  tip  deflections  knowing  the  strip  chart  strain 
gauge  readings . 
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EXPERIMENTAL  NATURAL  FREQUENCIES 
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LIST  OF  SYMBOLS 

A  Amplitude  of  oscillation 

ag  Two  dimensional  lift  curve  slope  =  2ir 

b  Semichord  =  c/2 

C(k)  Theodorsen  function 
Lift  coefficient 
c  Chord 

Flexibility  influence  coefficients 

Flexural  modulus  for  an  anisotropic  plate 

Et  Longitudinal  modulus 

Li 

Em  Transverse  modulus 

T 

e  Distance  between  elastic  axis  and  aerodynamic  center 

Ffc  Concentrated  test  load 

Gt  _  Shear  modulus 

I4X 

g  Structural  damping  coefficient 

i  /=T 

k  Reduced  frequency  =  wb/U 

k. .  Stiffness  influence  coefficient 
iD 

L_  Lift  force  at  elastic  axis 

E 

l  Length  of  plate 

Aerodynamic  Moment 
Mfc  Concentrated  test  moment 

m  Mass  per  area 


Distributed  lateral  load 


Generalized  force 
In-plane,  on-axis  lamina  modulus 
Dynamic  pressure 
Generalized  displacement 

Kinetic  energy 
Time 

Plate  thickness 
Velocity 

Divergence  velocity 
Flutter  velocity 

Strain  energy 

External  work 

Lateral  deflection 

Lateral  deflection  of  elastic  axis 

Tip  lateral  deflection 
P 

Complex  eigenvalue  =  (1  +  ig)/w2 

Angle  of  attack 
Initial  angle  of  attack 
Elastic  twist 
,  Tip  angle  of  attack 

Mode  shape  (two  dimensional) 

Ply  angle 


Coupling  factor 
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v__  Major  Poisson's  ratio 
vTL  Minor  Poisson's  ratio 
Air  density 

4> ^  Mode  shape  (one  dimensional) 

a)  Frequency  of  oscillation 
u)  Torsional  natural  frequency 


